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FOREWORD 


This  report  was  prepared  by  Dr.  Harold  H.  Sogia  and  Mr.  Kamal-Eldin 
Hassan  on  Contract  AF  33( 6l6 ) —98 •  They  were  assisted  by  Mr.  Frederick 
Salzberg.  The  work  was  performed  at  the  Heat  Transfer  Laboratory,  Depart¬ 
ment  of  Mechanical  Engineering,  Illinois  Institute  of  Technology,  Chicago, 
Illinois.  The  work  was  administered  by  Dr.  Max  Jakob,  Director,  and,  since 
his  demise  on  January  4»  1955*  by  Mr.  Stotbe  P.  Kezios,  Acting  Director  of 
the  Heat  Transfer  Laboratory.  The  contract  was  issued  from  the  Aeronautical 
Research  Laboratory,  Wright  Air  Development  Center  -  Mr.  Erich  Soehngen  was 
the  task  scientist  for  the  laboratory.  This  work  is  further  identified  as 
Task  70141,  "Transient  Heat  Transfer  in  Regenerators,"  under  Project  3066, 
"Gas  Turbine  Technology." 


ABSTRACT 


This  report  deals  with  the  thermal  design  of  rotary  type  regenera¬ 
tive  heat  exchangers  for  aircraft  gas  turbine  power  plants*  Hausen's 
regenerator  theory  is  developed  in  a  form  directly  applicable  to  the 
rotary  machine*  and  deviations  from  soma  of  the  underlying  assumptions 
are  examined.  A  numerical  method  to  calculate  the  performance  of  un¬ 
balanced  regenerators  is  included*  but  this  calculation  is  not  required 
for  application  of  the  theory.  Pressure  drop  and  heat  transfer  data  for 
flams  trap  and  wire  screen  matrices  are  assembled*  The  effect  of  the 
regenerator  performance  on  the  thermodynamic  cycle  is  discussed*  Tinally* 
a  sample  calculation  is  presented* 
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NOMENCLATURE 


Note  1:  The  symbol  F  is  used  for  the  unit  Fahrenheit  degree  (tenperature 
difference ) .  One  Rankine  degree  is  equal  to  one  Fahrenheit  degree; 
hence,  the  symbol  F  is  also  used  for  a  degree  on  the  absolute  tem¬ 
perature  scale.  The  symbol  *F  is  used  for  the  tenperature,  that 
is,  the  place  on  the  Fahrenheit  tenperature  scale. 

Note  2t  The  prime  (’),  the  double  prime  (") ,  and  the  triple  prime  (m)  are 
used  with  synfcols  C,  m,  q,  and  S  to  denote  that  the  quantity  refers 
to  a  unit  of  length,  a  unit  of  area,  or  a  unit  of  volume,  respectively. 


Symbol 

A 


C 

C,C1,C2,C3 

°P 

D 

D 

e 

f 

f 


G 

g 

g 

h 


h1>*«  *,hj^ 
K 
K 
k 
L 
l 


m 

m 


N 

N 


Quantity 

facial  area  or  cross  sectional  area  of  flow 
in  duct 

heat  capacity 

constants 

specific  heat 

wire  diameter 

equivalent  diameter  (see  Eq.  3-12) 

friction  factor  (see  Eq.  3-16) 

reduced  tenperature  of  the  solid  at  inlet  or 
outlet  plane  (see  Eq.  2-3h) 

mass  velocity 

acceleration  constant 

a  function 

coefficient  of  heat  transfer 
enthalpy 

fraction  of  circle  subtended  by  ducts 

kernel  functions  (see  Eq.  2-39  and  -I4I) 

thermal  conductivity 

matrix  length 

drum  length 

rate  of  flow 

index 

number  of  divisions  in  nurorical  calculation 
speed  of  rotation 
Nusselt  number,  hD^k 


Unit 


B/F 


B/lb  F 
in.,  ft 
ft 


It/ sec  ft2 
ft/ sec2 

B/hr  ft2  F 
B/lb 


B/hr  ft  F 
ft 
ft 

lb/hr 


sec 


-1 


rl 


Jr 


f 
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Symbol 


N. 


Pr 


N, 


Re 


N, 


St 

n 


n 

P 

<2a 

q 

R 


r 

r 

P 

S 

T 

t 

U 

u 

V 

v 

W 

x 


Y,Y 


y>y 


Quantity 

Prandtl  number,  ^  c^/k 
Reynolds  number,  OD^a  • p) 

Stanton  number,  a  •  h/Gc 

P 

exponent 

pitch 

pressure 

heat  added  during  the  thermodynamic  cycle 
rate  of  heat  transfer 
gas  constant 
radius 

pressure  ratio 

area  of  heat  transfer  surface 
absolute  temperature 
temperature  of  solid 
utilization  factor  (see  Eq.  2-20) 
fluid  velocity  in  x-direction 
volume,  relative  volume 
matrix  velocity  in  y-direction 
work  of  cycle 

distance  from  entrance  in  direction  of 
fluid  motion 

duct  width  (see  Fig.  2-2) 

distance  solid  particle  moves  in  stream 
(see  Fig.  2-2) 


Unit 


in."1,  ft 
lb/ft2 
B/lb 
B/hr 
ft/F 
ft 


*R 

*F 

ft/hr 

ft3 

ft/hr 

B/lb 

ft 

ft 

ft 


P 

Y 

A 

e 

\c 

^Rec 

^Heg 


ratio  of  free  flow  area  to  facial  area 

void  fraction 

specific  weight 

increment  of  £ 

dummy  variable  of  integration 

reduced  time  (see  Eq.  2-17) 

compressor  efficiency 

recuperator  effectiveness 

regenerator  effectiveness 
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Symbol 

Quantity 

Unit 

it 

turbine  efficiency 

— 

Q 

tenperature  of  gas 

*F 

It, 

ratio  of  specific  heats 

— 

A 

reduced  length  (see  Eq.  2-18) 

— 

/* 

dynamic  viscosity 

lb/hr  ft, 
lb/ sec  ft 

4 

reduced  distance  (see  Eq.  2-10) 

— 

IT 

reduced  period 

— 

P 

density 

lb/ft3 

T 

time 

hr 

<t> 

conduction  parameter  (see  Eq.  2-77) 

«— 

0 

variable  angle 

— 

10 

speed  of  rotation 

radians/ sec 

Subscripts 

A 

state  point  at  air-side  regenerator  outlet  in  thermodynamic 
cycle 

a 

air  side  of  regenerator 

app 

apparent 

B 

state  point  at  gae-side  regenerator  outlet  in  thermodynamic 
cycle 

b 

bulk  mean 

g 

gas  side  of  regenerator 

i 

inside 

m 

mean 

o 

entrance}  outside 

opt 

optimum 

r 

reduced}  radial 

s 

solid 

xiv 


i 


TRADC  TR  5^-13 


! 


Chapter  It  INTRODUCTION 


1-1  Purpose 

The  idea  of  regeneration  to  increase  the  thermal  efficiency  of  a 
thermodynamic  cycle  is  well  known.  Because  some  advantages  of  space 
and  weight  requirements  have  appeared  in  favor  of  the  rotating  regenera¬ 
tive  heat  exchanger,  many  investigators  during  the  past  decade  have 
examined  the  possibilities  of  its  application  to  gas  turbine  power 
plants  of  aircraft.  The  purpose  of  this  report  is  to  bring  the  present 
theoretical  and  experimental  knowledge  to  a  form  which  may  be  useful  to 
designers  interested  in  that  application. 

1-2  Scope 

The  report  deals  mainly  with  thermal  design  of  the  regenerator. 

The  influence  of  the  regenerator  performance  on  the  over-all  thermo¬ 
dynamic  cycle  is  considered.  Heat  transfer  and  pressure  drop  data  are 
presented.  Ifechanical  design,  power  plant  arrangement,  and  the  still 
inconpletely  solved  problem  of  satisfactorily  sealing  the  fluid  streams 
are  mentioned  only  insofar  as  they  may  have  a  bearing  on  the  thermal  de¬ 
sign  or  performance  calculations. 

1-3  Application  of  the  Thermodynamic  Regenerator 

In  the  next  few  sections  some  characteristics  of  thermodynamic  re¬ 
generation  are  briefly  reviewed  by  outlining  the  operation  of  a  gas  tur¬ 
bine  cycle  with  and  without  regeneration  and  by  conparing  the  performances 
of  the  two  cases. 

1-3.1  Thermodynamic  Cycle  of  the  Simple  Gas  Turbine  Power  Plant 

The  ideal  thermodynamic  cycle  of  a  sinple,  open,  gas  turbine 
power  plant  is  shown  diagranmatically  in  the  tenperature-entropy  plane 
of  Fig.  1-1.  The  air  at  the  atmospheric  condition  is  conpressed 

isentropically  to  state  p2,Tg.  111611  heat  is  add®d  at  constant  pressure 
by  burning  fuel  in  the  air  until  the  gas  attains  the  tenperature  T^. 

The  air  is  then  expanded  isentropically  in  a  turbine  to  atmospheric  pres- 
sur®  P^  =  Px* 

Manuscript  released  by  the  author  in  Juno  1955  for  publication  as  a  WOO 
Technical  Report* 
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1-3.2  Thermodynamic  Gas  Turbine  Cycle  with  Regeneration 

A  heat  exchanger  is  now  introduced  to  transfer  some  heat  from 
the  exhaust  gases  (Point  U)  to  the  high-pressure  air  flowing  from  the  com¬ 
pressor  to  the  combustion  chamber.  This  reduces  the  amount  of  fuel  needed 
to  bring  the  air  to  Point  3.  A  schematic  flow  diagram  of  a  regenerative 
gas  turbine  power  plant  is  shown  in  Pig.  1-2.  The  numbers  and  letters  re¬ 
fer  to  the  correspondingly  labeled  points  in  Pig.  1-1. 

As  can  be  seen  in  Fig.  1-1,  the  amount  of  heat  transferred  in 
the  exchanger  depends  upon  the  difference  between  the  temperatures  T^ 
and  Tg.  In  the  case  of  high  pressure  ratios  (Pg/p^),  it  can  happen  that 
Tli  ~^Z  so  that  a  heat  exchanger  would  decrease  the  thermal  efficiency. 
For  this  reason  regeneration  mist  not  be  introduced  arbitrarily. 

Farther,  with  the  addition  of  the  heat  exchanger  to  the  power 
plant  some  inherent  disadvantages  must  be  considered.  The  size  and  weight 
of  the  plant  increase,  the  pressure  losses  increase,  and  some  compressed 
air  leaks  to  the  exhaust  without  passing  through  the  turbine.  The  heat 
exchanger  must  do  more  than  just  compensate  for  these  deleterious  effects 
if  its  addition  is  to  be  worth  while. 

1-3.3  Evaluation  of  the  Heat  Exchanger  Effectiveness  in  the  Thermo¬ 
dynamic  Cycle 

The  ideal  amount  of  thermal  energy  that  could  be  transferred 
per  unit  time  from  the  exhaust  gases  to  the  relatively  cold  air  is 
maCp  ft(T^  -  Tg) ,  subscript  a  referring  to  the  cold  air  side  of  the  heat 
exchanger,  m  being  the  rate  of  flow  in  lb/hr,  and  cp  an  average  spe¬ 
cific  heat  in  B/lb  F.  For  practical  reasons,  the  actual  rate  of  heat 
transfer  is  m^c^  a(T^  “  Tg) ,  where  T^  is  the  actual  temperature  of 
the  cold  air  leaving  the  heat  exchanger  (see  Fig.  1-1  and  -2).  The  ratio 
of  the  actual  to  the  ideal  rate  is  called  the  heat  exchanger  effective¬ 
ness.  Accordingly, 

T  -  T 

^ex  =  T^-Tg  (1-1> 
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FIG.  1-2  REGENERATIVE  GAS  TURBINE  POWER  PLANT 


Schmitt  (Ref.  26)  calculated  the  influence  of  regeneration  on 
the  thermal  efficiency  of  a  gas  turbine  power  plant  operating  as  follows: 


Air  inlet  tenperature,  t^  =  -  58*F 
Maximum  cycle  tenperature,  t^  =  1832*F 
Heat  exchanger  effectiveness, 
Conpressor  efficiency, ^ 

Turbine  efficiency,^- 


or  Tj  =  l£2*R 
or  T3  =  22?2*R 

v =  0,7 

«]c  =  0,85 
It  =  °-82 


The  results  of  the  calculations  are  shown  graphically  in  Fig. 
1-3.  Obviously,  regeneration  is  most  practical  at  low  pressure  ratios, 
as  has  already  been  mentioned.  If  the  heat  exchanger  effectiveness  wuld 
be  increased,  the  maximum  thermal  efficiency  of  the  cycle  would  be  in¬ 
creased  and  translated  to  a  lower  value  of  the  pressure  ratio,  as  is  indi¬ 
cated  by  the  trend  of  the  curves  in  Fig.  1-3. 


One  way  to  increase  the  heat  exchanger  effectiveness  is  to  in¬ 
crease  its  size.  However,  there  is  an  optimum  value  of  n  beyond 
which  a  gain  in  cycle  efficiency  cannot  be  made  without  excessive  increases 
of  size  and  weight  of  the  power  plant.  This  optimum  depends  upon  the  heat 
exchanger. 


1-Lt  Some  Comparisons  of  Recuperative  and  Regenerative  Heat  Exchangers 

In  this  report  the  conventional  heat  exchanger  of  the  shell-and-tube 
type  is  called  a  "recipe  rat  or"  and  the  rotating  or  Ljungstroem  type  is 
called  a  "regenerator".  The  term,  "regenerator",  is  used  here  apart  from 
its  meaning  in  the  science  of  thermodynamics.  Where  ambiguity  may  arise 
the  longer  expression,  "regenerative  heat  exchanger",  is  employed. 

1-U.l  Limitation  of  the  Recuperator 

It  has  been  shown  that  the  smallest  recuperators  to  meet  a  spe¬ 
cified  performance  are  obtained  by  using  large  nurfcers  of  tubes  having 
small  diameter  (see,  for  exanple,  Ref.  28).  One  of  the  main  reasons  is 
that  in  fine  passages  the  flow  becomes  laminar,  in  which  case  the  required 

^The  cospressor  and  turbine  efficiencies  are  defined  in  Chapter  I4. 
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FIG.  1-3  EFFICIENCIES  OF  SIMPLE  AND 
REGENERATIVE  GAS  TURBINE 
CYCLES 
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core  size  diminishes  such  more  rapidly  with  the  Reynolds  number  than  in 
the  case  of  turbulent  flow.  A  practical  limit  of  the  tube  size  is  set 
by  constructional  difficulties  in  separating  the  two  fluids  inside  the 
recuperator  core  and  in  providing  headers  or  tube,  sheets;  also,  operational 
difficulties  arise  on  account  of  fouling  and  cleaning. 

1-U.2  Regenerator  Core  or  Matrix 

In  the  regenerator  core  the  constructional  difficulties  can  be 
avoided,  ffell  defined  passages  are  unnecessary,  for  the  two  fluids  are 
allowed  to  occupy  the  same  space  at  different  times.  The  heat  is  exchanged 
primarily  by  virtue  of  the  matrix  capacity  for  heat  storage  rather  than 
by  its  nature  to  conduct  heat.  Hence,  the  size  of  the  passages  nay  be  con¬ 
siderably  reduced,  and  a  corpact  and  relatively  light  weight  core,  or  matrix, 
can  be  obtained.  Since  seals  within  the  matrix  are  unnecessary,  any  hetero¬ 
geneously  packed,  permeable  solid  can  be  employed.  Thus,  the  matrix  may 
be  constructed  from  layers  of  wire  gauze  or  screens,  odd  shapes  of  braided 
wires  and  ceramic  fibers,  as  well  as  various  types  and  shapes  of  crimped 
ribbons  and  tubular  packings,  There  is  some  evidence  (Ref.  1 5)  indicating 
that  no  appreciable  amounts  of  soot  build  up  inside  of  the  passages  of 
these  matrices,  and  so  the  above  mentioned  operational  difficulties  may 
be  avoided. 

1-5  Regenerator  Types 

The  matrix  of  the  regenerator  is  exposed  to  the  hot  fluid  for  some 
time  during  which  its  temperature  rises.  Then  it  is  exposed  to  the  cold 
fluid  which  takes  away  the  amount  of  heat  that  had  been  stored  during  the 
heating  period.  This  process  repeats  itself  in  a  cyclic  way.  In  appli¬ 
cations  of  the  steel  and  liquefaction  industries,  the  regenerators  are 
built  with  two  matrices  which  are  alternately  heated  and  cooled  by  switch¬ 
ing  the  fluids  with  valves.  The  important  theories  of  regenerators  have 
been  developed  with  those  arrangements  in  view. 

It  seems  that  the  method  most  suitable  for  producing  this  cyclic 
change  in  aircraft  power  plants  is  to  rotate  the  matrix  so  that  any  part 
of  the  matrix  is  carried  alternately  through  the  hot  and  cold  fluid  streams. 


1ADC  TR  55-13 


7 


This  method  is  similar  to  that  of  the  rotative  type  air  preheaters  used 
in  modern  steam  power  plants.  In  this  report  one  of  the  main  regenerator 
theories  is  recast  so  that  it  is  directly  applicable  to  the  rotating  re¬ 
generator  j  also,  other  theories  are  discussed  from  that  standpoint. 

Rotating  regenerators  nay  be  classified  according  to  Aether  the 
fluids  flow  axially  or  radially. 


1-5.1  Axial  or  Disk  Regenerators 

The  matrix  of  the  axial  regenerator  shown  in  Fig.  1-1*  is  built 
in  the  shape  of  a  disk.  The  disk  is  mounted  coaxially  in  a  round  cylinder 
divided  into  two  ducts.  The  two  fluids  flow  in  opposite  sense,  the  hot 
fluid  flowing  in  one  duct  and  the  cold  fluid  in  the  other.  The  matrix  ro¬ 
tates  about  its  axis  and  would  usually  be  turned  by  a  positive  external 
drive  requiring  a  very  small  amount  of  power.  Thus,  any  point  in  the 
nfitrix  passes  alternately  through  the  hot  and  cold  streams,  carrying  the 
heat  from  the  one  to  the  other. 


1-5.2  Radial  or  Drum  Regenerators 

Figure  1-5  represents  the  drum  type  regenerator.  The  matrix  is 
built  in  a  cylindrical  form  through  which  the  fluids  flow  radially.  This 
drum  is  mounted  coaxially  with  two  round  ducts.  The  inner  duct  and  the 
annular  space  are  divided  in  the  axial  direction  into  two  parts,  one  for 
the  cold  fluid  and  the  other  for  the  hot  fluid.  As  in  the  previous  case, 
the  fluids  are  in  counterflow  through  the  matrix. 

1-6  Leakage 

One  of  the  practical  problems  arising  in  the  mechanical  design  of  a 
regenerator  for  application  to  the  gas  turbine  power  plant  is  the  prevention 
of  leakage,  particularly  from  the  high  pressure  air  to  the  low  pressure 
gas.  The  amounts  of  leakage  play  an  inport  ant  role  in  the  thermal  design 
and  thermodynamic  performance.  Four  types  of  leakage  may  possibly  occur. 
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FIG.  I  -  4  DISK  TYPE  REGENERATOR 


FIG. I* 5  DRUM  TYPE  REGENERATOR 
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1-6.1  Air-to-Gas  Leakage  through  Clearance  Spaces 

Leakage  from  the  air  side  to  the  gas  side  occurs  through  gaps 
between  the  moving  matrix  and  the  stationary  walls  separating  the  two 
fluids.  This  type  of  leakage  affects  the  power  plant  efficiency  adverse¬ 
ly,  because  the  compressor  must  handle  air  at  a  higher  rate  than  the  tur¬ 
bine.  The  leakage  can  be  considerably  reduced  by  the  inside  seals  shown 
in  Fig.  1— U  and  -5. 

1-6.2  Air-to-Alr  or  Gas-to-Gas  Leakage 

Some  air  or  gas  nay  bypass  the  matrix  through  the  clearance 
space  between  the  rotating  matrix  and  the  duct  walls,  because  a  pressure 
drop  exists  across  the  matrix.  Hence,  parts  of  either  fluid  may  not  pass 
through  the  matrix,  and  the  effectiveness  of  the  regenerator  diminishes.. 
This  type  of  leakage  is  reduced  by  means  of  the  outside  seals  indicated 
in  Fig.  1— U  and  -5. 

1-6.3  Air-to-Gas  Leakage  through  Matrix 

In  matrices  other  than  the  tubular,  or  flame  trap  type,  the 
air  could  leak. through  the  matrix  to  the  exhaust  side.  This  is  prevented 
by  dividing  the  matrix  into  sectoral  c (apartments  and  using  seals  wide 
enough  to  cover  at  least  one  compartment  at  a  time. 

1-6. U  Carry-Over  Leakage 

Finally,  the  so-called  "carry-over"  leakage  or  "let  down"  oc¬ 
curs.  The  fluid  trapped  in  the  individual  conpartments  of  the  matrix  are 
carried  from  one  stream  to  the  other  by  virtue  of  the  matrix  rotation. 
Because  the  air  has  a  high  density,  the  carry-over  is  greater  from  the  air 
to  the  gas  side  than  from  the  gas  to  the  air  side  of  the  regenerator.  The 
carry-over  is  unavoidable;  however,  it  can  be  kept  below  a  reasonable 
limit  with  slow  speeds  of  rotation. 


1-7  Other  Design  Considerations 

In  addition  to  the  leakage  which  is  a  problem  unique  to  the  regenera¬ 
tive  heat  exchanger,  the  pressure  drop,  common  to  all  types  of  heat  ex- 
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changers,  oust  be  considered.  The  effects  of  the  pressure  drops  on  both 
the  air  and  gas  sides  is  to  decrease  the  expansion  ratio  of  the  turbine, 
thus  reducing  the  thermal  efficiency  of  the  power  plant.  To  be  worth 
while,  the  regenerator  mist  of  course  augment  the  net  output  considerably 
more  than  would  be  needed  just  to  balance  the  decrease  due  to  pressure 
drop. 

The  effect  of  size  and  weight  should  be  considered  from  the  viewpoint 
of  the  aircraft  performance  and  purpose.  As  already  mentioned,  because 
of  its  small  passages  and  finely  divided  structure,  the  regenerator  appears 
to  have  an  advantage  of  being  lighter  than  a  recuperator  meeting  the  same 
performance  requirements. 
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Chapter  2 :  THEORY  07  THE  ROTATING  REGENERATIVE  HEAT  EXCHANGER 

* 

2-1  Assumptions 

The  assumptions  of  Hausen’s  (Ref.  10  and  12)  theory  are  discussed  in 
the  next  sections.  Then  his  basic  differential  equations  are  developed 
from  the  standpoint  of  ttieir  application  to  the  rotating  regenerator. 


2-1.1  Homogeneity 

It  is  assumed  that  both  the  solid  and  the  fluids  are  continuous, 
homogeneous  media.  This  assumption  is  essential  for  the  formulation  of 
the  differential  equations.  In  practice  the  matrix  may  be  considered 
homogeneous  if  the  material  is  so  finely  divided  that  the  si2es  of  the 
voids  and  individual  solid  elements  are  very  small  compared  to  the  over¬ 
all  volume  of  the  matrix. 


2-1 . 2  Steady-State  Operation 

It  is  assumed  that  the  matrix  has  been  rotating  long  enough  so 
that  the  temperature  of  any  point  in  the'  matrix  solid  repeats  itself  peri¬ 
odically.  Thus,  any  transients'  of  starting  and  stopping  are  not  considered. 

At  each  point  in  space  through  which  the  matrix  passes  the  tem¬ 
perature  is  independent  of  the  time.  This  means  that  as  an  element  in 
the  matrix  reaches  a  particular  position  in  space,  it  attains  a  certain 
temperature  affixed  to  that  position. 


2-1.3  Uniform  Temperatures  of  the  fluids  at  Entrances  to  the  Matrix 
The  temperature  of  the  fluids  at  the  entrances  to  the  matrix  is 
assumed  uniform.  Of  course,  at  the  exits  their  temperatures  vary  from 
point  to  point  as  will  be  determined  by  the  theory. 

2-1. b  Convection  without  Conduction 

Heat  is  transferred  only  by  convection  at  the  interface  of  the 
matrix  and  the  fluids.  An  infinitesimal  elemental  subdivision  of  the 
matrix  material  interchanges  heat  only  with  the  contacting  fluid  and  not 
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with  neighboring  subdivisions j  within  each  infinitesimal  subdivision  of  the 
solid  the  conductivity  is  infinite  in  the  direction  normal  to  the  direction 
of  fluid  flow  and  zero  in  the  direction  parallel  to  the  fluid  flow.  In 
other  words,  no  heat  is  conducted  through  any  finite  part  of  the  solid, 
while  the  capacity  of  this  part  for  heat  storage  is  not  impaired.  Since 
the  matrix  is  usually  made  of  thin  sheets  or  wires  of  heat  resisting  ma¬ 
terials  which  have  relatively  low  conductivity,  this  assumption  is  justi¬ 
fiable.  Some  effects  of  finite  conductivity  are  discussed  in  later  sections 

As  for  the  fluids,  heat  is  convected  by  them  and  is  conducted 
within  them  only  in  regions  on  the  interface  of  each  subdivision.  This 
assumption  is  justified  by  the  low  conductivity  of  the  fluids. 


2-1.5  Constant  Coefficients  of  Heat  Transfer 

The  coefficient  of  heat  transfer  h  is  assumed  constant  in  either 
stream.  This  means  that  influences  of  transients  in  the  boundary  layers 
and  of  variations  in  the  fluid  properties  are  neglected.  This  assumption 
may  be  justified  by  the  facts  that  the  transients  are  of  relatively  brief 
duration  and  that  the  thermal  properties  of  the  fluids  change  relatively 
little  with  tenperature.  In  practice,  the  fluid  properties  are  evaluated 
at  certain  mean  temperatures  (see  Section  5-6). 

2-1.6  Constant  Density 

The  thermal  expansion  or  contraction  of  the  fluid  during  heating 
or  cooling  changes  the'  velocity  of  the  fluid  as  it  passes  through  the 
matrix.  The  influence  of  these  variations  are  omitted  from  the  theory  in 
order  to  simplify  the  mathematics.  Calculations  by  Saunders  and  Snoleniec 
(Ref.  25)  indicate  that  this  assumption  is  justifiable. 

2-2  Derivation  of  the  Differential  Equations 

Suppose  that  the  matrix  has  been  stretched  out  so  that  its  motion  is 
rectilinear.  Then  the  matrix  is  said  to  be  flat  or  plane.  It  is  evident 
that  the  findings  regarding  the  flat  matrix  are  directly  applicable  to  the 
disk  type  regenerator j  also,  it  is  shown  in  Appendix  A  that  after  minor 
modification  the  results  are  applicable  to  the  drum  type  regenerator. 


WADC  TR  55-13 


13 


With  reference  to  Fig,  2-1,  the  matrix  has  the  length  L  in  the  main 
direction  of  fluid  flow.  The  velocity  of  the  matrix  is  v,  the  direction 
being  normal  to  that  of  the  fluid  stream.  The  fluid  flows  with  a  mass 
velocity  of  m"  pounds  per  hour  per  square  foot  of  facial  or  frontal  area. 
A  heat  balance  is  made  on  an  element  of  volume  in  space.  The  dimensions 
of  the  element  in  the  plane  of  motion  are  5x  by  8y$  its  depth,  normal 
to  the  plane  of  motion,  is  one  foot. 

Under  steady-state  conditions  (see  Section  2-1.2),  thermal  energy  is 
carried  in  the  x-direction  by  the  fluid  and  in  the  y-direction  by  the  solid 
also  some  heat  is  carried  over  in  the  y-direction  by  the  fluid  trapped  in 
the  interstices  of  the  matrix  element. 

2-2.1  Thermal  Energy  Conveyed  by  the  fluid 

The  energy  transported  per  unit  time  by  the  fluid  into  the  ele¬ 
ment  of  space  is 


a  =  m"  (By  •  1)  •  c  •  6  (2-1) 

jC  p 

where  0  =  0(x,y)  is  the  fluid  temper ature }  c^  •  0  is  the  enthalpy  of 
the  fluid  referred  to  zero  enthalpy  at  0*F. 


2-2*2  Thermal  Energy  Conveyed  by  the  Solid  Material 

Similarly  as  in  Section  2-2.1,  the  thermal  energy  transported 
per  unit  time  by  the  solid  material  into  the  element  of  space  is 


^,1  ■  T  •  (8x  '  1}  ’  fs.app  •  cp,s  •  *  <2-« 

where  subscript  s  refers  to  the  solid,  t  =  t(x,y)  is  the  solid  tempera¬ 
ture,  and  .c  •  t  is  the  enthalpy  of  the  material  referred  to  zero  en- 
p,s 

thalpy  at  O  F,  Subscript  app  means  that  the  apparent  density  should  be 
enployed.  Subscript  1  indicates  that  only  one  part  of  the  heat  conveyed 
in  the  y-direction  is  considered. 


11*  . 
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FtG;2-l  NOMENCLATURE  FOR  THE  PLANE 
MATRIX  REGENERATOR 
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2-2.3  Thermal  Energy  in  the  Carry-Over 

A  second  part  of  the  heat  conveyed  in  the  y-direction  is  the  therm¬ 
al  energy  of  the  fluid  carried  over  in  unit  time  by  the  matrix* 


V.8  "  VitX  •  1}  ?*PP 


(2-3) 


2-2.1*  Balance  of  the  Transported  Thermal  Energy 

Applying  the  First  Law  of  Thermodynamics  to  the  elemental  volume 
of  space  and  observing  that  no  heat  is  stored, 

^2E  5x  +  ^V*1  *  Vi.2!  5y  as  0  (2-1*) 

dx  dy 

Substituting  for  q^,  q^.  and  q^.  ^  ^**om  Eq.  2-1,  -2,  and  -3  into  Eq.  2-1* 
and  simplifying, 

a"  c  —  +  C"»  v  —  +  p  fl  vc  —  =  0  (2-5) 

p  dx  dy  '  '  p  dy 

where  ^  is  the  actual  fluid  density,  p  is  the  void  fraction  (the  ratio 
of  the  volume  of  voids  to  the  volume  of  matrix),  and 


C*  sp  .  c 

“s,app  p,s 


(2-6) 


2-2.5  Heat  Exchange  in  the  Elemental  Volume  of  Space 

As  the  relatively  hot  matrix  "flows"  into  the  space  element  it 
loses  heat  to  the  flowing  fluid  by  convection.  The  rate  of  change  of 
the  enthalpy  of  the  matrix  material  is  equal  to  the  rate  of  heat  trans¬ 
ferred  by  convection  in  the  same  space  element.  Accordingly, 

^  (v  •  5x  •  C*  t)  5y  =  -  hSw  (t  -  ©)  5x  5y  (2-7) 

The  minus  sign  appears  because  t  decreases  in  the  positive  y-direction. 
The  same  equation  would  be  obtained  if  the  matrix  were  considered  to  be 
cold  and  gaining  heat  from  the  fluid.  Simplifying, 


* 
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(2-8) 


at  hs"1 

"Sy  “  vCw 


(0  -  t) 


Elininating  dt/dy  from  Eq.  2-5  and  -8  and  placing  m"  =  upp, 


90 

Sc 


+ 


V 

u 


90 

9y 


(t  -  0) 


(2-9) 


It  nay  be  noted  that  u  is  the  effective  axial  velocity  of  the  fluid,  that 
is,  the  distance  a  fluid  particle  moves  per  unit  time  in  the  x-direction 
through  the  natrix  interstices.  Equations  2-8  and  -9  are  the  basic  dif¬ 
ferential  equations  governing  the  tenperatures  of  the  solid  and  fluids  in 
the  flat  regenerator.  In  the  next  section  they  are  brought  to  a  more  con¬ 
venient  form. 


2-2.6  The  Differential  Equations  in  Dimensionless  Space  Coordinates 
Introducing 

5  =  “ffF  x  (2-3.0) 

m"c_ 


n  -  *5  $  -  f>  (2-u) 

and  substituting  for  x  and  y  in  terms  of  5  and  ,  Eq.  2-8  and  -9 
become 

(2-12) 

(2-13) 

where  t  and  0  are  now  regarded  as  functions  of  £  and  ^  • 

Equations  2-12  and  -13  describe  the  heat  exchange  in  one  of  the  re¬ 
generator  streams.  Therefore,  the  problem  of  the  heat  exchange  in  the 
regenerator  is  solved  using  two  sets  of  these  equations.  In  both  streams, 


2£  =  e-t 

dr, 

*  =  t-e 

95, 
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x  is  positive  in  the  direction  of  the  floe  and  y  is  positive  in  the  di¬ 
rection  of  the  matrix  notion* 


2-3  Boundary  Conditions 

Only  counterflow  regenerators  will  be  considered.^  The  two  streams 
will  be  called  the  "cold  air"  and  the  "hot  gas",  and  symbols  denoting 


quantities  on  corresponding  sides  of  the  regenerator  will  be  distinguished 
by  subscripts  a  and  g,  respectively.  Thus,  a  plane  normal  to  the  flow 


direction  of  the  cold  air  is  at  a  distance  x  from  the  cold-air  entrance. 

Cl 

The  same  plane  would  be  at  a  distance  x^  from  the  hot-gas  entrance,  so 
that 


xg  =  L  -  \  (2-lli) 

In  order  to  formulate  the  boundary  conditions  the  "actual"  regenerator 
will  be  replaced  by  an  "ideal"  regenerator.  It  will  be  convenient  to  con¬ 
sider  the  matrix  of  fine  tubular  construction;  however,  substantially  the 
same  analysis  would  apply  to  gauze  type  matrices. 


2-3.1  Difficulty  of  Formulating  Boundary  Conditions  for  the  Actual 


As  the  solid  particles  in  a  plane  normal  to  the  direction  of  the 
matrix  motion  enter  one  of  the  fluid  streams,  a  point  at  x^  ^  0  is  not 
immediately  affected  because  the  fluid  does  not  reach  that  point  until 
the  time  x^/u  has  elapsed;  meanwhile  the  point  moves  to  y^  srv  • 

Also,  at  the  farther  side  of  the  fluid  stream,  at  x  =  0,  there  is  a  dis¬ 
tance  in  the  y-direction,  namely,  v  «  (L/u),  such  that  fluid  particles 
entering  cannot  reach  all  parts  of  the  passages  before  they  are  blocked 
or  sealed.  Summarizing,  in  the  actual  regenerator  the  duration  of  time 
a  point  in  the  matrix  is  affected  by  a  given  stream  depends  upon  its  dis¬ 
tance  x  from  the  entrance  of  the  fluid.  A  point  of  the  matrix  at 

I  1 

x  =  0  spends  the  full  heating  or  cooling  period  in  the  fluid  stream 
whereas  a  point  at  x  =  L  stays  a  shorter  time.  These  variable  exposure 

^"Hausen  (Ref.  10)  also  considered  parallel-flow  regenerators.  He 
showed  that  their  maximum  effectiveness  is  limited  to  a  value  of  0.5. 
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tines  lead  to  mathematical  boundary  conditions  which  are  difficult  to 
handle.  Therefore,  an  ideal  regenerator  is  defined  which  closely  approxi¬ 
mates  the  actual  regenerator  for  all  practical  purposes  and  simplifies 
the  mathematics. 


2-3.2  Ideal  Regenerator 

An  ideal  flat  regenerator  is  shown  diagrammatically  in  Fig,  2-2. 
The  walls  dividing  the  two  fluid  streams  are  imagined  to  be  permeable, 
in  such  a  way  that  they  allow  the  matrix  to  pass  through  them  but  keep 
the  fluids  out.  Further,  the. walls  of  each  channel  are  slanted  at  the 
angle  tan-^  (v/u)  from  the  main  direction  of  the  fluid  flow.  Hence,  the 
fluid  particles  move  parallel  to  the  walls.  As  a  consequence  of  this 
arrangement,  any  point  in  the  matrix  enters  the  stream  at  the  moment  the 
fluid  arrives;  the  point  is  immediately  affected  by  the  first  particle 
of  fluid  entering  at  the  near  side  of  the  stream.  At  the  far  side,  the 
fluid  particles  leave  the  matrix,  without  being  replaced  by  other  fluid 
particles.  Thus,  in  the  ideal  regenerator,  carry-over  from  one  stream 
to  the  other  is  avoided. 


With  reference  to  Fig.  2-2, 


and  the  distance  between  the  permeable  walls  is 


(2-15) 


Y*  =  T  -  l  x  (2-16) 

Hence,  the  distances  y  are  measured  from  an  axis  -which  is  oblique  to 
the  direction  of  flow,  and  the  distance  Y  is  uniform,  so  that  all 
points  of  the  matrix  in  the  ideal  regenerator  have  the  same  length  of 
effective  path  in  a  given  stream. 

It  may  be  noted  that  upon  eliminating  y  from  Eq.  2-11  and  -15, 


hS"'  « 

n  =  7 


(2-17) 
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FIG. 2-2  IDEAL  REGENERATOR 
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and  that  the  differential  equations  in  dimensionless  space  coordinates, 
namely,  Eq.  2-12  and  -13,  are  unaltered. 


2-3*3  Dimensionless  Lengths 

To  formulate  the  boundary  conditions  of  the  ideal  regenerator, 
the  extreme  values  of  S.  and  trill  be  needed.  The  value  of  £  for 
x  =  L  trill  be  denoted  by  A  and  trill  be  called  the  reduced  length.  Thus, 

A  S  L  (2-18) 

P 

# 

The  value  of  for  y  =  Y  trill  be  denoted  by  It  and  tail  be  called 

the  reduced  period.^-  Thus, 


TI 


=  hS™  * 
~  vCm  1 


(2-19) 


It  is  of  some  interest  that  the  ratio 


I) 


m"c 


i  =  ^-T 


(2-20) 


■which  is  called  the  utilization  factor,  is  independent  of  h  and  S"'. 


It  should  be  noted  that  values  of  A.  and  H  on  the  air  side 
of  the  regenerator  are  generally  not  equal  to  corresponding  values  on  the 
gas  side.  Since  x  is  measured  from  the  plane  at  which  the  fluid  enters 
the  matrix  and  since  s/A  =  x/L,  the  quantities  £  and  t  corre- 
sponding  to  the  same  plane  are  related  by  the  equation. 


(2-21) 


The  reason  for  this  bizarre  name  for  a  distance  is  that  the  theory 
was  originally  developed  for  non-rotary  regenerators-  which  operate  in  a 
quasi-stationary  way,  and  the  term  ■period"  as  well  as  the  symbol  TI 
are  retained  to  avoid  any  confusion  in  case  the  reader  consults  the 
literature . 
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Further, 


A  m? 


_£1a 


m"  c 

g  P»g 


!k 

h. 


and 


•* 

Tl  Y 

— £  =  -£ 

Tl.  v* 


•* 


Eliminating  the  ratio  h g/hft  from  Eq.  2-22  and  -23, 


ae -  . 

A." 


m"  c  Y 

a  p,a  a 

* 

*■  c  Y 

g  P»g  g 


(2-22) 


(2-23) 


(2-2U) 


Moreover, 


For  many  practical  purposes  it  nay  be  assumed  that  c  _  =  c  , 

p,a  p,g 

in  the  application  of  the  theory  to  an  aircraft  gas  turbine  power  plant, 
if  the  leakage  from  one  stream  to  the  other  is  very  small  cospared  with 

f  # 

the  total  rate  of  flow  it  may  also  be  assumed  that  m"  Y  »  m"  Y  ,  by  coi 

a  a  g  g»  J 

tinuity.  Under  these  circumstances  the  ratio  in  the  parentheses  of  Eq, 
2-2li  takes  the  value  1,  and 


=  Jfi  or  V  =  n  ■  (2-25) 

A.  1.  a  * 

This  is  a  useful  relationship  in  dealing  with  regenerators  which  have  un¬ 
equal  reduced  lengths  and  reduced  periods,  that  is,  so-called  unbalanced 
regenerators. 


2-3.1;  Boundary  Conditions  of  the.  Ideal  Regenerator 

In  accordance  with  the  assumption  in  Section  2-1.3  that  the  tem¬ 
peratures  of  the  fluids  at  the  entrance  planes  are  uniform, 

V0'  “  ®a,0  <2-26> 

and 


V°*  V  “  ®g»o 
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where  both  0.  _  and  0  „  are  constants. 
a,o  g,o 

Since  nothing  happens  to  the  ideal  matrix  as  it  passes  through 

the  permeable  wall,  the  temperature  at  a  point  in  the  solid  leaving  one 

stream  is  equal  to  the  teuperature  at  that  point  entering  the  second 

stream.  Accordingly,  considering  that  t  =  t  (x  ,y  )  and  t  =t  (x  ,y  ), 

8  £  g  6  a  a  a  a 

and  focusing  our  attention  on  a  plane  of  the  matrix  which  leaves  the  cold 
air  and  on  another  which  enters  the  hot  gas,  we  may  write 


W»  -  VL  -  y<> 


(2-28) 


Introducing  the  reduced  lengths  end  reduced  periods  and  considering  now 
that  tg  =  tg(  and  ta  =  t&(  vja)>  we  may  transform  Eq.  2-28 
to 


e„>°>  =  » 


(2-29) 


Similarly,  focusing  our  attention  on  planes  which  leave  the  hot  gas  and 
enter  the  cold  air. 


v  v>  =  *«  [*  -  -£]•  y)  (2-30> 

The  problem  of  the  ideal  regenerator  is  now  mathematically 
fonnulated.  A  solution  is  required  to  satisfy  two  sets  of  the  linear 
partial  differential  Eq.  2-12  and  -13 >  one  set  for  each  fluid,  and  to 
satisfy  the  boundary  conditions,  namely,  Eq.  2-26,  -27*  -29,  and  -30. 

2-3.5  Conparison  of  the  Ideal  and  Actual  Regenerators 

The  difference  lies  in  the  slanted,  permeable  wall,  which  is 
impossible  to  construct  in  the  actual  regenerator.  This  wall  has  two 
features: 

first,  it  has  the  ability  to  separate  the  fluid  from  the  solid, 
allowing  the  solid  to  pass  through  it.  The  effect,  as  already  mentioned, 
is  to  eliminate  the  carry-over  losses  discussed  in  Section  1-6.  ij.  In 


WADC  TR  55-13 


23 


practice  this  type  of  leakage  may  be  kept  low  by  employing  low  speeds  of 
rotation,  say,  30  to  1$  revolutions  per  minute.  This  places  a  limit  on 
the  matrix  velocity  v. 

Second,  the  permeable  wall  is  slanted  at  the  angle  tan“^  (v/u) 
with  respect  to  the  direction  of  fluid  flow,  while  the  wall  of  the  actual 
regenerator  is  straight.  As  (v/u)  -»  0,  the  shape  of  the  ideal  regenera¬ 
tor  approaches,  the  shape  of  the  actual  one.  Since  practical  values  of 
v/u  would  usually  be  quite  small,  the  approximation  is  good.  It  may  be 
noted  that  the  same  approximation  would  be  made  by  omitting  the  second 
term  in  the  parentheses  of  Eq.  2-11,  for  this  term  arose  by  the  admittance 
of  a  carry-over  of  thermal  energy  from  one  elemental  volume  of  space  to 
another  (Section  2-2.3).^ 


2-3.6  Single-Blow  Problem 

An  allied  problem,  which  is  of  a  sinpler  nature  than  the  one 
treated  above,  deals  with  the  "single-blow".  Its  inportance  here  lies 
in  the  fact  that  its  solution  may  be  employed  to  solve  the  regenerator 
problem.  Also,  it  has  been  enployed  in  several  experimental  studies  to 
determine  the  coefficients  of  heat  transfer.  See,  for  exanple.  Ref.  12, 
25,  and  26. 

In  the  single-blow,  a  stationary  matrix  having  an  arbitrary 
initial  tenperature  distribution  is  suddenly  subjected  to  a  continuous 
blow  of  fluid  entering  with  constant  tenperature.  Mathematically,  the 
problem  is  to  solve  one  set  of  the  basic  differential  equations,  namely, 
Eq.  2-12  and  -13,  with  one  boundary  of  the  form  of  Eq.  2-26  or  -27  and 
the  condition  that 


t(5,0)  =  f(e ) 


(2-31) 


where  f {  £  )  is  a  prescribed  function.  The  solution  is  conprised  of  ex¬ 
pressions  for  t(  5  ,  n  )  and  ©(  £  ,  rj ) ,  Solutions  are  presented  in 


The  next  several  sections  deal  with  the  mathematical  solution  of  the 
regenerator  problem.  The  reader  who  is  interested  mainly  in  the  results 
and  their  application  may  proceed  without  loss  of  continuity  to  Section  2-6 . 
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Ref*  2,  22,  and  27;  the  last  two  references  deal  with  the  case  that 
f C  4  )  =  constant. 

2— It  Solutions  of  the  Regenerator  Problem 

Hausen  (Ref.  10)  developed  a  purely  analytical  solution.  The  main 
results,  presented  in  Section  2-8,  are  directly  applicable  to  the  balanced 
regenerator,  that  is,  a  regenerator  in  which  =  Ag  and  Ua  = 

Nusselt  (Ref.  23),  working  independently,  brought  the  mathematical  problem 
to  a  form  requiring  the  solution  of  a  set  of  integral  equations.  Numeri¬ 
cal  methods,  one  of  which  is  based  on  Nusselt's  treatment  and  is  presented 
in  Section  2-5,  have  since  shown  that  the  results  for  the  balanced  re¬ 
generator  may  be  applied  to  predict  the  over-all  performance  of  an  un¬ 
balanced  regenerator,  provided  a  proper  averaging  procedure  is  employed. 
Although  quite  time  consuming,  the  numerical  methods  have  the  advantages 
of  being  applicable  to  both  balanced  and  unbalanced  regenerators  and  of 
providing  the  transient  histories  and  the  steady-state  temperature  distri¬ 
butions  in  both  the  matrix  and  the  fluids.  Even  when  the  tenperature 
distributions  of  a  balanced  regenerator  are  required,  a  numerical  solution 
appears  preferable  to  direct  substitution  into  the  mathematical  series 
comprising  the  results  of  Hausen's  analytical  solution. 

Hausen  (Ref.  11)  devised  the  first  numerical  method,  the  so-called 
heat-pole  method.  He  divided  the  matrix  into  a  number  of  strips  or  heat 
poles  along  the  x-axis.  Since  no  conduction  occurred,  he  could  treat 
each  pole  as  though  it  were  an  individual  bed  in  the  single-blow  theory, 
for  which  analytical  solutions  of  a  closed  form  were  already  available. 

Saunders  and  Smoleniec  (Ref.  25)  solved .the  single-blow  problem 
numerically  by  converting  the  differential  equations  to  difference  equa¬ 
tions  and  integrating  in  a  stepwise  manner.  Starting  with  a  regenerator 
matrix  of  uniform  tenperature,  allowing  it  to  cycle  successively  through 
the  hot  and  cold  fluids,  and  observing  the  point  in  the  calculations 
where  the  tenperatures  repeated  themselves,  they  determined  the  steady- 
state  distributions  and,  consequently,  the  regenerator  performance.  At 
the  same  time  they  could  predict  the  number  of  revolutions  the  matrix 
turned  before  coming  to  steady-state  conditions.  In  any  practical  case 
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where 


Ka<  6 )  5  “ 


_  «L(avvs.-‘ 


(2-39) 


In  the  same  way,  Eq.  2-35  becomes 

i 

l-f.  A»[x-i!]  ■•"‘Nv] 


■fvv’h'*’] 


where. 


W*5"-  -v 


«i  (  2iV  Tt„<  6.  - 


{g-  * 


(2-1*0) 


(2-ia) 


Equations  2-38  and  -1*0  are  two  simultaneous  integral  equations 
giving  the  temperature  distributions  in  the  matrix  planes  at  the  entry  to 
the  two  fluids,  namely,  f  and  f  .  It  will  be  seen  that  when  these 
distributions  are  known,  the  effectiveness  can  readily  be  found.  To 
carry  out  the  calculations,  the  integrals  are  to  be  evaluated  by  means  of 
Simpson's  rule  or  its  equivalent,  as  shown  below. 


2-5.2  Review  of  Simpson’s  Rule  for  Integration 
2*ne*,!*  "  '  "P  _ 

To  evaluate  a  definite  integral  /  g(e)  de  by  Simpson's  rule, 

**  8l 

the  interval  (a,b)  on  the  6 -axis  is  divided  into  an  even  number  of  parts, 
each  of  width  A  .  Approximating  the  function  g(.e  )  with  parabolic  seg¬ 
ments  through  every  group  of  three  successive  points,  one  may  write  the 

2 

following  equation  with  an  error  of  the  order  C\. 
rb 

/  g(  &  )  d*  =  -y  £g(a)  +  ^g(a  +  A  )  +  2g(a  +  2  A)  +  i*g(a  3  A) 


+  2g(a  -f  1*A)  +  ...  +  2g(b  -  2 A) 


+  l*g(b  -  A  )  + 


g(b)J 


(2-1*2 )] 


The  parentheses  in  this  equation  and  in  all  remaining  equations  of 
this  section  enclose  the  functional  arguments. 
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In  particular)  if  the  number  of  intervals  is  2, 


j: 


g(e)  de 


a  +  b 


+ 


(2-1*3) 


2 -5.3  Application  of  Sjnpson's  Rule  to  the  Solution  of  the  Integral 

Equations 

In  applying  Sinpson’s  rule  to  solve  Eq.  2 -38  and  -1*0,  each  of 
the  reduced  lengths  A  and  A  is  divided  into  N  equal  parts  of 

a  g 

length  and  Ag»  respectively,  as  shown  in  Fig.  2-3.  It  is  unneces¬ 

sary  that  N  be  an  even  number.  Regarding  Eq.  2-38,  consider  that 
£  =  m  •  A  where  m  =  0,  1,  2,  . . .,  N.  Then 

CL  6L 


fg(n  •  Ag  0  -  jp>  =  ra<*  •  Aa)  •  e 


-It. 


I 


m  •  A 


Ka(m  •  Aa“  & )  •  fa(fe  )  de 


(2-1*1*) 


With  the  understanding  that  a  symbol  of  the  type  f  will  henceforth 

8>n 

replace  the  symbol  of  the  type  f  (n  •  A-)*  Eq.  2-1*1*  takes  the  sinpler 

©  o 

form. 


f  =  f 

g,N-  m  a,m 


-V 

e  + 


m 

/ 


Ka(m 


^a" 


e) 


fa< 


e)  de 


(2-1*5) 


This  equation  is  typical  of  N  +  1  equations  which  will  be  written  by 
placing  m  =  0,1,2,...,  and  N.  Equation  2-1*0  could  be  written  in  a 
similar,  typical  manner  to  represent  another  set  of  N  +  1  equations. 
Employing  SLnp son’s  rule  or  its  equivalent  to  evaluate  the  integrals,  one 
obtains  a  total  of  2N  +  2  linear  algebraic  equations  in  which  the  un¬ 
knowns  are  f  and  f  for  m=  0,1, 2,... ,N. 

g,m  a, a  7  ’ 


Since  Sinpson's  rule  can  be  used  only  if  m  is  an  even  number, 
equivalent  expressions  are  needed  for  the  cases  in  which  m  is  odd.  In 
the  next  four  sections,  the  integrals  are  developed  according  to  whether 
m  is  even,  a  =  1,  m  =  3,  and  a  is  an  odd  nunber  greater  than  3. 
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2-5.1*  Integration  for  a  Even 

In  this  case  Sinpson’s  ruLe  is  directly  applicable  and  the  inte¬ 
gral  in  Eq.  2-1*5  becomes 


/ 


m*  A 


K  (m  •  A.  -  e)  •  £„(&)  d€r 

cl  Cl  ft 


Ka(m  •  Aa-0)  •  fa(0)  +  l*Ka(m.  Aft-  •  f&(  Aa) 


+  2Ka(m  •  Aa-  2  Aa)  *  fa(2  Aa)  +  •  •  • 

+  2Ka(m*  Aa-[a-2]  Aa)  *  fa([®“2].  Aa) 
+  l*Ka(m-  Aa-[m-l]  Aa)  •  fa ([»-!]*  Aft) 


+  Ka(m  •  Aa  -  a  ‘  Aa)  *  fa(“  •  Aft) 


(2-1*6) 


The  symbol  K  may  replace  the  lengthier  notation  K  (n  •  A.)>  a  simi- 

81  y  n  ft  ft 

lar  change  of  notation  having  already  been  made  in  regard  to  f.  Then 

m  •  A 
*-*a 


/. 


Ka(m  •  Aa-  e)  •  fa(  «  ) 


T 1 


K  •  f  _  +  IjK  ,  •  f  ,  +  2K  •  f  . 

a,m  a,0  ^a,m^l  a,l  a,m-2  a, 2 


+  1*K  _  ,'f  ,+  .„+2K  « 


+  1*K 


.  ,  •  f  .  +  K  .  •  f 
a,l  a,m-l  a,0  a,m 


(2-1*7) 


In  the  special  case  that  m  =  0,  the  integral  is  zero. 


2-5.5  Integration  for  m  =  1 

In  this  case  Iliffe  assumed  that  the  tenperature  distribution 
can  be  Represented  in  the  range  0  ^  s-  ^3  A  by  a  polynomial  of  the 
third  order: 


V  e  )  “  *0  +  *lfe  +  a2  e  2  +  a^  e  ■ 


(2-1*8) 
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Substituting  e  =  0,  Aa,  2Aa,  *nd  3Aft>  the  value8  of  *0*  *1*  a2* 

and  a,  can  be  determined  in  terms  of  f  f  - ,  f  0,  and  f.  Then 
i  &*x  a*£ 

enploying  Eq.  2— 1*8  to  determine  fa^2  s  fa(  Aa/2), 


f 

a 


1 

A/2  =  13 


+  l$f 


a,l 


(2—1*9) 


Thus,  the  interval  (0,  A  )  has  been  subdivided  into  two  equal  parts  so 

a 

that  Eq.  2-1*3  is  applicable.  The  result  is 


/ 


K.(Aa  -  6)  •  f,U)  d. 


T"  1  Ka,l  '  fa,0  +  £  Ka,l/2 

•  r#.  rt+i5f ,  -f  0 +f  ,i 

L  a>°  aA  a*2  a»3j 


+  "a.O  •  f.,l 


(2-50) 


2-5.6  Integration  for  m  =  3 « 

Iliffe  assumed  that  the  entire  integrand  could  be  represented 
by  another  third-order  polynomial: 


Ka(3Aa  -  e  )  •  fa(e)  =  aq  +  A^  +  Ag  e2  +  k^  fe3 


(2-51) 


Again  substituting  fe  =  0,  A  ,  2  A  ,  and  3  A  ,  the  values  of  A  ,  A_ , 

a  a  a  Ox 

k0,  and  A-  can  be  determined  in  terms  of  K  •  f  K  •  f  , 

£  3  a^O  a,_L 

K  _  •  f  and  K  •  f  After  replacing  the  coefficients  in  Eq.  2-51 

with  their  newly  derived  values  and  integrating  with  respect  to  e- 

f rom  0  to  3  A  » 

a 


3A„ 


/  Ka(3 Aa”  e)  •  fa(e)  de 
‘'o 

■  1  A.  [Ka,3  •  f«,0  +  3K.,2  •  f.,l+  3K.,1  •  fa,2  +K.,0'  fa,3]  (2‘52) 
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2-5.7  Integration  for  m  Odd  and  Greater  than  3 
In  this  case  Iliffe  eirployed  the  relationship 

^a  r3  *  Aa  Aa 

I  g(  e )  de  a  /  g(  &  )  de  +  /  g(  e  )  de 

^0  Jo 


(2-53) 


Thus,  he  could  use  the  results  from  the  case  that  m  =  3  in 
conjunction  with  the  result  from  the  case  that  m  is  even.  With  due  re¬ 
gard  for  the  actual  meaning  of  the  index  m  in  the  present  case,  for 
exanple,  observing  that  m-3  is  even,  the  final  result  is 

rm- A. 

/  Ka(m  •  Aa-  &  )  •  fa(  e)  d e 

Jo 

=  ^a  |  F  Ka,m  *  ^a,0+F  ^a,m-l  *  ^a,l  Ka,m-  2  *  ^a,2 
+  [f+t]  Ka,m-3  *  fa,3+I  Ka,m-  h  'fa,i*+f  Ka,m-5  *  fa,5  + 


+  4  K 


—  K  _«f  „+-=K  i  •  f  ,  +4  K 

3  a, 2  a,m- 2  3  a,l  a,n-l  3  a,0  a,mf 


(2-5W 


2-5.8  Final  Form  of  the  Simultaneous  Linear  Algebraic  Equations 

Applying  the  results  of  Sections  2-5. U  through  2-5.7  successively 

in  accordance-  with  Eq.  2-U5  for  m  =  0,1,2,...,N,  one  obtains  the  N  +  1 

equations  shown  in  Fig.  2-h;  they  are  referred  to  as  Eq.  2-55.  As  already 

suggested,  the  other  set  of  N  +  1  equations  my  be  obtained  by  replacing 

f  ,T  ,  A  >  K  >  f  >  and  11  by  1  -  f  „  >  A  f  K  ,1-f  , 
g ,.N  -  nr  ^a*  a, nr  a,m  iVa.  J  a,N-nr  ug’  g,nr  g^ 

and  ft  ,  respectively.  These  equations  are  shown  in  Fig.  2-5  and  are 

o 

referred  to  as  Eq.  2-56. 


2-5.9  General  Remarks  regarding  the  Calculation  Procedure 

The  sets  of  Eq.  2-55  and  -56  are  2N  +  2  linear  algebraic  equa¬ 
tions  in  f  and  f  for  m=  0,1, 2,... ,N.  The  solution,  therefore. 

g,m  a,m  ’  * 

represents  the  tenperature  distributions  in  the  matrix  at  the  planes  of 
entry  to  the  gas  and  the  air. 
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FIG.  2-4  EQUATIONS  2-55 


+  e” 
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FIG.  2-5  EQUATIONS  2-56 


GENERALLY  FOR  ANY  EVEN  VALUE 
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FIG.  2-5  (CONT.)  EQUATIONS  2-56 


The  functions  K  defined  by  Eq.  2-39  and  —III  always  take  posi¬ 
tive  values  and  can  be  evaluated  by  means  of  tables  of  Bessel  functions. 

_  2 

For  values  of  2yit{  £  -  &}  greater  than  I,  the  following  formula 

can  be  used  to  determine  J^(ix)* 


Jx(ix) 


ie 


y  5  ir  x 


1  - 


0.375  0.1171875  0.10253906  0.1M1956 


x3 


jr 


0.2775761 

“7" 


(2-57) 


Once  the  distributions  of  f  and  f  are  established  the  re- 

a  g 

generator  effectiveness  can  be  calculated  as  shown  in  Section  2-7*  Also, 
the  temperature  at  any  point  in  the  matrix  can  be  determined  by  means  of 
Eq.  2-33  and  -35.  Further,  the  temperature  distributions  of  the  fluids 
can  be  calculated  by  means  of  Eq.  2-12. 


2-6  Numerical  Calculation  of  Balanced  Regenerators 

In  the  case  that  A  =  A_  =  A  and  TC  =  Tl  =  TC  a  symmetry 
occurs  in  the  two  sets  of  algebraic  equations  of  Section  2-5.8  which 
makes  them  identical,  namely, 


f  w  =  1  -  f 
g,N-m  a,m 


(2-58) 


In  this  case  it  is  sufficient  to  solve  just  one  set  of  equations  and  to 
employ  Eq.  2-58. 


2-7  Regenerator  Effectiveness  obtained  from  Numerical  Calculations 

In  this  section  the  effectiveness  is  calculated  from  the  viewpoint 
of  the  heat  capacity  of  the  matrix.  The  matrix  enters  the  cold  stream 
with  a  certain  enthalpy  and  leaves  it  to  enter  the  hot  stream  with  a  lower 
enthalpy.  In  steady  conditions,  the  difference  between  these  two  enthalpies 


^See,  for  example,  E.  Jahnke  and  F.  Erode,  Tables  of  Functions  with 
Formulae  and  Curves,  Dover  Publications,  New  Yorlc7"T^3715TT57^  r~ 

2Ibid . ,  pp.  137-138. 
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is  the  actual  amount  of  heat  energy  transferred  from  the  matrix  to  the  cold 
stream;  an  equal  amount  of  energy  is  picked  up  by  the  matrix  as  it  tra¬ 
verses  the  hot  stream. 


The  ideal  amount  of  heat  energy  transferred  is  the  quantity  that  the 
cold  air  would  pick  up  if  its  temperature  had  been  raised  to  that  of  the 
incoming  hot  gas.  The  ratio  of  the  actual  total  amount  to  the  ideal  total 
amount  is  the  regenerator  effectiveness.  The  following  calculations  are 
made  with  the  above  procedure  in  mind. 


For  a  unit  width  of  matrix,  the  rate  at  which  heat  is  carried  by  the 
matrix  as  it  moves  into  the  cold  stream  is 


Q 


g 


-4 -I 

,o  a,o 


vC"  fa(  £a)  ^ 


mHc 

=  vC«»  .  P.il  / 

u  h  S"  / 
a  JQ 


fa<  d*. 


(2-59) 


Similarly,  the  rate  at  which  heat  is  carried  by  the  matrix  as  it  moves 
out  of  the  cold  air  is 

»  Agr 

(2-60) 


q’ 


m"c 


6  -  0 
g,o  a,o 


=  vC"'  •  r-£- 


g 


p  I 


v  «*>  d  «« 


The  ideal  rate  of  heat  transfer  would  be  obtained  if  the  tenperature 

anj 


of  the  cold  air  would  change  from  ©a  Q  to  ®g  D*  Accordingly, 


0  -  Q 

g,o  Wa,o 


=  m"  c 


a  p,a 

Hence,  the  regenerator  effectiveness  is 


(2-61) 


^Reg  q7" 

,  r\  n  m"  c  Y*  rAg 

■* :l  vv*«, 

The  integrals  in  Eq.  2-62  are  conveniently  evaluated  in  terms  of  the 

numerical  values  of  f  and  f  obtained  from  the  solution  outlined  in 

6 
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Section  2-5*  Applying  Sinp son's  rule  for  the  case  that  the  number  of 
intervals  N  is  even, 

rA 


f  tu  )  dl  =  Ff0  +  ^  +  2f2  +  iif3 


+  2f  ^  +  •  •  • 


+  2fN-2  +  lrfN-l+fN 


(2-63) 


And  if  N  is  odd,  the  technique  of  Section  2-5*7  may  be  employed  -with 
the  result  that 

=y[|f0-f  *i+f  *htk*uS 


+  k f ^  +  2f j  +  •••  +  2f _  g 


^N-l  +  fJ 


(2-61i) 


2-8  Results  of  the  Theory  of  the  Balanced  Regenerators 

The  results  of  Hausen's  analytical  solution  or  of  the  corresponding 
numerical  calculations  are  most  conveniently  represented  graphically. 

Figures  2-6  and  -7  give  the  effectiveness  of  balanced  regenerators 
for  values  of  A  and  It  in  the  range  of  practical  interest.  The 
curves  were  drawn  from  data  given  by  Johnson  (Ref.  15).  In  Fig.  2-8  curves 
of  constant  utilization  factor  are  presented.  The  graphs  illustrate  the 
following  properties  of  regenerators: 

(1)  For  a  given  effectiveness  there  is  a  minimum  reduced  length  A 
below  which  that  effectiveness  cannot  be  attained. 

(2)  An  increase  in  the  reduced  length  A  (or  actual  length  L)  is 

more  effective  in  the  case  of  a  regenerator  with  a  small  value  of  A 

than  one  with  a  large  value  of  A  . 

(3)  For  a  given  value  of  A  »  the  effectiveness  decreases  as  It 

increases.  This  effect  is  quite  large  for  small  values  of  A  • 
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REGENERATOR  EFFECTIVENESS 


FIG. 2-6  EFFECTIVENESS  OF  BALANCED 

REGENERATORS 

Uo 
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REDUCED  LENGTH  A 

FIG. 2-7  EFFECTIVENESS  OF  BALANCED 
REGENERATORS 

41 
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2-9  Unbalanced  Regenerators 

The  results  of  Section  2-8  can  be  used  to  obtain  the  effectiveness 
of  unbalanced  regenerators,  at  least  in  the  range  of  practical  applica¬ 
tion.  A  method  of  procedure  suggested  by  Saunders  and  Smoleniec  (Ref.  25), 
namely,  to  enploy  a  system  of  arithmetic  averages,  has  been  found  to  yield 
results  -which  deviate  from  the  exact  values.  Two  methods  of  procedure, 
based  on  harmonic  averages,  have  been  found  to  give  virtually  equal  values 
and  better  accuracy  when  conpared  with  results  determined  by  the  much  more 
lengthy  procedure  described  in  Sections  2-5  and  -7.  In  the  next  two 
sections  it  is  supposed  that  the  parameters  A  and  It  on  both  sides 
of  the  regenerator  are  known. 


2-9.1  Use  of  Harmonic  Means  of  Dimensionless  Lengths 

First,  the  harmonic  rneetns  of  A  and  A  and  of  1  and 

a  g  a 

Ti.  are  determined : 

g 


2 

Am 


(2-65) 


Second,  employing  charts  for  the  balanced  regenerator  find 

the  calculated  values  of  A  and  H_.  The  effectiveness 

m 


1n»g,m  using 

m'  ^Begjm 

practically  equal  to  the  effectiveness  of  the  given  unbalanced  regenerator. 


2-9 . 2  Use  of  Harmonic  Wean  of  Effectiveness 

Hausen  suggested  the  following  procedure.  Determine  the  effect¬ 
iveness  r^g  a  of  a  balanced  regenerator  with  parameters  A&  and 
and  the  effectiveness  ^ Reg  g  R  balanced  regenerator  with  parameters 

A  and  Tt  .  Then  the  effectiveness  of  the  given  unbalanced  regenerator 
&  & 

is  practically  the  harmonic  mean.  Thus, 


1 

^Regja 


1  • 

^Reg,g 


(2-66) 


2-10  Effect  of  Curvature  in  the  Drum-Type  Matrix 

In  Appendix  A  it  is  shown  that  a  sinple  transformation  allows  all  re¬ 
sults  of  the  calculations  for  a  flat  matrix'  to  be  employed  in  dealing  with 


h3 


WADC  TR  55-13 


a  drum-type  matrix.  The  major  changes  occur  in  the  transformations  of 
Section  2-2.6,  where  t  and  v|  are  defined.  The  new  values  of  S. 
and  ^  are 


£  = 


hSw 

m*c 

P 


2ir  K 


2 


(2-67) 


and 


hS*  0 
1  =  cw  L 


2-trKpp 


(2-68) 


In  these  equations  m'  is  the  mass  rate  of  fluid  flow  per  unit  length  of 
the  regenerator  axis;  K  is  the  fraction  of  the  circle  subtended  by  the 
fluid  in  question;  r  is  a  general  radial  distance;  r^  is  the  inner 
radius  of  the  drum;  0  is  the  angle  measured  from  the  plane  of  entrance; 
and  w  is  the  angular  speed.  These  coordinates  allow  for  the  changes 
of  tenperature  distributions  caused  by  the  curvature  of  the  matrix. 


The  corresponding  new  values  of  A  and  IT  are: 


and 


A  = 


hS™ 

m"c 
m  p 


L 


(2-69) 


j  -  hS-  .  2j[K 

11  ~  U3 


(2-70) 


where  L  is  the  thickness  of  the  drum  and  is  the  mass  velocity 
of  the  fluid  at  the  mean  radius  of  the  drum. 

2-11  Simplified  Theories  for  Regenerators  with  Small  Reduced  Periods 

In  aircraft  applications,  the  reduced  periods  for  the  regenerators 
are  rather  short,  of  the  order  of  magnitude  of  5.  The  short  reduced 
periods  can  be  due  either  to  high  matrix  speed  v  or  to  the  small  ratio 
of  the  water  equivalents  of  the  flowing  fluids  to  that  of  the  "flowing" 
matrix. 
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For  infinitely  short  periods,  i.e.,  for  H  =  K  =0.  Nusselt 
(Ref.  22)  showed  that  the ' regenerator  would  have  the  same  thermal  charac¬ 
teristics  as  a  recuperator.  Under  this  condition,  the  matrix  tenperature 
depends  only  on  x  (or  d  ),  is  independent  of  y  (or  ),  and  is  the 

same  in  the  two  streams.  Also,  the  temperatures  of  the  two  fluids  are 
functions  of  x  only.  For  this  problem,  Nusselt  considered  both  the 
balanced  and  unbalanced  cases. 

Iliffe  (Ref.  13)  further  showed  that  for  unbalanced  regenerators  with 

infinitely  short  periods,  if  m"  c  Y  =  mM  c  Y  ,  then  the  effective- 

a  P>a  a  g  P>g  g 

ness  would  be  that  of  an  unbalanced  recuperator,  namely, 

W  =  r—x 

Aa  Ag 

For  balanced  regenerators  this  reduces  to  the  balanced  recuperator 
effectiveness, 


^Rec  A +  2 


(2-72) 


According  to  Tipler  (Ref.  29),  the  problem  of  the  high-speed  re¬ 
generator  with  linear  temperature  distributions  of  the  matrix  was  solved 
by  Lubbock  and  Bowen  (Ref.  19).  The  general  assumptions  of  Hausen's 
theory  are  the  basis  of  the  solution.  Tipler  represents  the  results  as 
follows  s 


Aa 

*a 


’iReg 


=  coth 


*a 


+  coth 


g 


(2-73) 


Tipler  compared  values  of  thus  obtained  v/ith  those  from  Hausen's 

theory  for  balanced  regenerators.  He  found  that  for  II  =  0  the  agree¬ 
ment  is  good  and  that  for  other  values  the  effectiveness  derived  by  this 
method  is  less  than  Hausen's  values.  For  values  of  TT  up  to  10  the 
error  is  within  U  per  cent.  Up  to  It  =  5,  which  covers  a  large  range  Of 
practical  interest,  the  maximum  error  is  1  per  cent. 
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Corbitt  (see  discussion  of  Ref.  13)  presents  the  sinple  equation, 

Vg  =  'W  (1  ’  T5  (2'7U) 

where  ^Rec  va^ue  obtained  from  Eq.  2-71  or  -72.  This  relation 

appears  to  be  accurate  within  1  per  cent  in  the  range  of  ^Reg  ^rom  0*8 
to  0.98  and  of  U  from  0  to  0.6. 

2-12  Effect  of  Conductivity  on  Regenerator  Performance 

The  assunption  of  Section  2-1. Ij,  which  specifies  that  the  matrix  be 
made  of  a  certain  type  of  thermally  non-isotropic  solid,  cannot  be  exactly 
fulfilled  in  practice.  In  the  next  sections  the  influences  of  conduction 
in  the  matrix  are  considered. 


2-12 .1  Conductivity  in  the  x-Direction 

When  the  general  procedure  of  Section  2-2  is  followed  and  the 
conductivity  in  the  x-direction  is  admitted  to  the  problem,  Eq.  2-5  ap¬ 
pears  with  another  term: 

2 

m"c  jg  +  C"  v  ^  +  p  (3v  c~-  k(l  -  (3)  =  0  (2-75) 

P  dy  \  '  P  dy  r 

?here  k  is  the  effective  conductivity  of  the  matrix  in  the  x-direction. 
Also,  Eq.  2-7  would  be  replaced  by 

2 

vC»"  ~  =  -  hS"  (t  -  0)  +  k(l  -  G)  (2-76) 

?>y  r  ^2 

Equation  2-9  would,  therefore,  remain  unchanged  and  together  with  Eq.  2-7 6 
would  determine  the  rate  of  heat  transfer  in  the  regenerator.  It  appears 
that  this  general  problem  has  never  been  solved  analytically. 

However,  the  case  of  infinitely  large  conductivity  in  the  x- 
direction  was  solved  by  Nusselt  (Ref  .  22).  In  this  case  the  natrix  tem¬ 
perature  is  independent  of  x  and  the  maximum  effectiveness  is  0.5. 

This  theory  could  possibly  be  applied  to  very  short  matrices. 
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Also,  Hahnemann  (Ref.  8)  solved  the  special  problem  of  the  case 
that  TV  =  ft  =  0.  Here  the  heat  capacity  is  not  a  criterion  of  per** 

O  | 

formancej  it  may  take  any  value  without  affecting  values  of  It  or  A  . 
Hahnemann  found  that  to  represent  the  results  he  could  use  the  parameter, 

,  _  hS"L  ^ 

*  =  grpr  (2‘77) 

which  represents  the  ratio  between  the  heat  transfer  by  convection  and 
that  by  conduction.  It  may  be  observed  that  increases  as  k  de¬ 

creases,  so  that  for  a  matrix  with  zero  conduction  in  the  x-direction, 
as  assumed  by  Hausen,  <t>  =  co. 

The  essential  results  of  Hahnemann's  work  are  shown  in  Fig.  2-9.* 
Hausen's  matrix  is  represented  by  the  straight,  diagonal  line  for  4>  =  oo. 
It  is  observed  that  as  $  decreases  (i.e.,  k  increases),  the  regenerator 
effectiveness  decreases.  For  a  given  value  of  $  ,  there  is  an  optimum 
value  of  A  ,  at  which  the  regenerator  effectiveness  is  a  maximum. 

The  last  observation  may  be  explained  as  follows:  The  general 
tendency,  as  shown  in  Pig.  2-9,  is  that  r^-363  with  increasing  A  . 

However,  if  A  is  already  large,  then  the  rate  of  change  of  ^Rgg 
with  respect  to  A  is  rather  small.  Now,  suppose  that  A  is  already 
large  and  is  increased  by  diminishing  m",  say,  by  increasing  the  super¬ 
ficial  cross-sectional  area  of  flow.  This  change  would  have  no  influence 
on  <J>  .  However,  the  heat  conduction  per  unit  rate  of  flow  would  increase, 
and  at  a  certain  value  of  A  ,  namely,  the  optimum  value,  this  loss  due  to 
the  increase  of  conduction  offsets  the  gain  due  to  the  increase  of  A. 

As  A  is  increased  further,  the  conduction  losses  become  excessive  and 
diminishes. 

Another  parameter  used  in  Fig.  2-9  is  the  quantity, 

A/$  5  k(l  -  p)  /  (mw  0^  L).  It  may  be  considered  as  the  ratio  of  the 
rate  of  heat  conduction  along  the  natrix  to  the  rate  at  which  heat  is 
carried  by  the  fluid.  Cox  end  Stevens  (Ref.  5)  used  Hahnemann's  results 
for  balanced  regenerators  to  plot  Fig.  2-10,  which  covers  the  range  of 
interest  in  regenerator  design.  Again,  the  parameter  is  A/d>.  This 

*  This  figure  was  deleted  from  this  report.  It  is  available  in  the 
literature  cited  as  Reference  8. 
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NOTE 

This  figure  has  been  deleted 
from  this  technical  report,  but 
is  available  in  the  literature 
cited  as  Reference  &  of  this 
report# 
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FIG  2-10  EFFECT  OF  CONDUCTION  ALONG  THE  MATRIX 
ON  THE  EFFECTIVENESS  OF  THERMALLY 
BALANCED  REGENERATORS  WITH  0-0 


parameter  would  seldom  be  more  than  0.02,  which  in  the  usual  range  of 
operation  reduces  the  regenerator  effectiveness  by  less  than  0.015*  For 
most  design  purposes  the  effect  of  conductivity  in  the  matrix  could  be 
neglected.  It  is  advisable  to  check  this  effect  when  dealing  with  very 
short  matrices. 

2-12.2  Conductivity  in  the  z-Direction 

In  the  main  theory  the  matrix  was  assumed  homogeneous  (see 
Section  2-1.1).  This  implies  that  solid  particles,  as  well  as  voids, 
have  infinitesimal  volumes.  Under  such  conditions,  each  solid  particle 
would  be  affected  directly  by  the  flowing  fluids  while  neighboring  solid 
particles  would  not  affect  each  other.  A  real  matrix,  however,  consists 
of  solid  elements  of  finite  size,  and  some  solid  particles  do  not  come  in 
direct  contact  with  the  fluids  but  exchange  heat  by  conduction  through 
neighboring  particles.  For  exarple,  if  one  considers  a  matrix  built  of 
solid  strips  as  shown  in  Fig.  2-11,  then  the  particles  at  the  center  of 
&  strip  are  heated  and  cooled  by  virtue  of  the  heat  conduction  through 
particles  lying  between  the  center  and  the  surface.  If  those  strips  are 
of  infinitesimal  thickness,  then  each  solid  particle  would  exchange  heat 
directly  with  the  flowing  fluids  and  the  original  assumptions  would  still 
be  satisfied.  Also,  if  the  conduction  of  the  strip  in  the  z-direction  is 
infinitely  large  (while  that  in  the  x-  and  y-directions  are  zero),  then 
the  original  assumptions  are  also  satisfied,  because  there  would  be  no 
resistance  to  heat  flow  between  the  flowing  fluids  and  the  inner  solid 
particles. 

The  case  of  a  matrix  with  finite  conductivity  in  the  z-direction 
and  zero  conductivity  in  the  x-  and  y-directions  was  considered  by  Nus- 
selt  (Ref.  22).  He  gave  a  graphical  method  of  solution.  Ackermann 
(Ref.  1),  a  student  of  Nusselt,  undertook  an  analytical  solution  and  de¬ 
veloped  an  iterative  method  which  is  very  tedious. 

The  effect  of  finite  conductivity  in  the  z-direction  was  esti¬ 
mated  by  Tipler  (Ref.  29)  in  the  following  manner:  He  considered  the 
case  of  a  slab  whose  surface  temperatures  change  sinusoidally.  He  calcu¬ 
lated  the  heat  stored  in  the  slab  at  the  peak  surface  temperature  and 
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FIG. 2-11  REGENERATOR  WITH  MATRIX 

COMPOSED  OF  RIBBONS 

/ 


corpared  it  -with  the  heat  that  would  be  stored  in  the  slab  had  it  been  of 
infinite  conductivity.  He  found  that  the  ratio  of  the  two  amounts  would 
be  0.99  or  more  for  practical  regenerators.  Hence,  for  all  practical 
purposes,  the  effect  of  finite  conductivity  in  the  z-direction  can  be 
neglected. 


Nusselt  (Ref.  22)  also  gave  a  solution  for  the  case  of  a  strip 
matrix  with  infinitely  large  conductivity  in  the  x-direction,  zero  con¬ 
ductivity  in  the  y-direction,  and  finite  conductivity  in  the  z-direction. 
The  solution  was  given  as  expressions  of  t  and  0  in  Fourier  series 
in  terms  of  y  and  z. 

2-12.3  Conductivity  in  the  y-Direction 

In  the  main  theory  the  effective  conductivity  in  the  y-direction 
is  assumed  to  be  zero.  The  case  of  finite  conductivity  in  that  direction 
seems  never  to  have  been  considered  in  the  literature.  However,  the  case 
of  infinitely  large  effective  conductivity  in  the  y-direction  is  equivalent 
to  that  of  infinitely  high  matrix  speed  v,  which  has  been  discussed  in 
Section  2-11. 


2-13  Effect  of  Leakage  on  Regenerator  Performance 

The  theoretical  treatments  are  based  on  the  assunption  that  no  leakage 
occurs.  In  practice  the  leakage  diminishes  the  effectiveness  by  cooling 
the  hot  gas  before  its  heat  capacity  can  be  most  fully  utilized  in  the 
matrix. 

Harper  and  Rohsenow  (Ref.  9)  calculated  the  reduction  in  the  re¬ 
generator  effectiveness  assuming  that  half  the  leakage  occurs  at  either 
side  of  the  matrix.  Their  results  are  shown  in  Fig.  2-12  for  several 
values  of  fractional  leakage  One  may  employ  this  chart  to  obtain 

the  effectiveness  a  regenerator  with  leakage.  If  it  were 

assumed  that  all  leakage  occurs  at  the  downstream  air-side  of  the  matrix 
(i.e.,  leakage  into  the  gas  at  the  upstream  end  of  the  gas-side)  the  cor¬ 
rection  would  be  nearly  doubled,  with  further  reduction  of  the  effective¬ 
ness.  However,  it  is  more  probable  that  the  leakage  would  be  less  on  the 


YADC  TR  55-13 


52 


CORRECTION 


J0.06 


0.02 


POOL _ 1 _ 1 _ 1— 

— 

■  i 

0.20  0.40 

0. 

60  O.i 

:alculated  regenera 

>R  effectin 

12  EFFECT  OF  LE 

A 

KAGE  ON 

REGENERATOR 

E 

FFECTIVEI 

downstream  air-side  than  on  the  upstream  side.  Consequently,  an  effective¬ 
ness  calculated  by  means  of  the  chart  in  Fig.  2-12  would  be  slightly  low. 

For  the  practical  ?<vige  of  effectiveness  and  for  tolerable  values  of 
the  leakage,  the  correction  is  about  0.01,  which  is  not  serious.  However, 
the  effect  of  leakage  on  the  thermal  efficiency  of  the  entire  gas  turbine 
power  plant  is  a  serious  natter.  This  subject  receives  detailed  attention 
in  Chapter  U. 
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Chapter  3:  PROPERTIES  OF  REGENERATOR  MATRICES 
3-1  Void  and  Solid  Fractions 

Several  geometric  properties  and  some  of  the  mechanical  and  thermal 
properties  of  matrices  are  readily  expressed  in  terms  of  the  void  fraction, 
which  is  denoted  by  p.  The  co-value,  the  solid  fraction  (1  -  p),  of  a 
given  matrix  of  mass  mg  and  corresponding  total  volume  Is  the 

ratio  of  the  apparent  density  to  the  true  density  of  the  solid: 

(1  -  =  r-V  =  5-¥ee  (>1) 

*  tot  rs  fs 

It  is  reconmended  that  Eq.  3-1  be  enployed  as  the  operational  definition 
of  a  rather  sinple  laboratory  procedure  to  determine  p  for  matrices  of 
heterogeneous  construction  having  unknown  packability. 


3-1 .1  Flame  Trap  Matrices 

The  void  fraction  of  flame  trap  matrices  is  equal  to  the  ratio 
of  through  area  to  facial  area.  This  ratio  is  easily  c alculated  if  the 
passages  repeat  themselves  in  any  uniform  way. 


3-1 . 2  Wire  Screen  Ma trices 

Consider  screens  made  of  -wires  having  uniform  diameter  D  ex¬ 
pressed  in  feet.  The  number  of  wires  per  foot  in  one  direction  is  n^ 
and  in  the  other  n0j  the  screens  are  stacked  so  that  the  number  of  layers 
per  foot  of  length  is  n^j  then  neglecting  any  bends  of  the  wire,  the 
solid  fraction  is 

(1  -  p)  =  -J-  D2  +  n2)  n3  (3-2) 

If  the  screens  are  closely  packed,  n^  l/(2D)  and 

(1  -  (3)  =  -f-  0(1^  +  n2)  (3-3) 

Further,  if  n^  =  ng  =  n,  then 

(1  “  P)  “  T  D  *  n  (M) 
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3-2  Heat  Transfer  Surface  Area  per  Unit  Volume 

The  effective  heat  transfer  surface  per  unit  volume  of  matrix  is 


S»n 


(3-5) 


The  last  expression  allows  Sm  to  be  evaluated  in  terms  of  the  mass  of 
the  solid  material  used  to  build  the  matrix,  if  it  can  be  assumed  that 
the  contact  areas  are  negligible  conpared  with  the  areas  exposed  to  the 
fluids. 


For  exanple,  consider  a  matrix  made  of  any  metal  strips  of  uniform 
rectangular  cross  section;  let  a  be  the  thickness  and  b  the  breadth 
of  the  strips.  If  the  total  length  of  the  strip  is  \  ,  the  total  sur¬ 
face  area  is  2(a  +  b)  \  ,  and  its  mass  is  ab  •  \  .  ^)g.  Consequently, 
for  such  a  matrix  Eq.  3-5  yields 

S»  =  -(1  “#■  (1  +  S)  (3-6) 

a  b 


And  if  the  strip  is  very  thin,  so  that  a  <gib, 

S«  =  (3-7) 

CL 


In  a  similar  way,  if  wire  of  diameter  D  is  used  to  construct  the 
matrix,  Eq.  3-5  becomes 

S«  =  --CLl  (3-8) 

In  particular,  if  the  wires  are  in  the  form  of  closely  packed  screens  of 
square  mesh  (n  by  n  wires  per  foot). 


Sw  =  irn 


(3-9) 


3-3  Free  Area  of  Fluid  flow 

The  free  flow  area  is  needed  to  evaluate  the  velocity  through  the 
matrix  interstices,  which  is  of  interest  in  determining  the  pressure  drop 
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and  heat  transfer  characteristics  of  the  matrix.  In  the  following  dis¬ 
cussion  it  will  be  convenient  to  denote  the  ratio  of  the  free  flow  area 
to  the  facial  area  by  the  symbol  <x. 

In  a  flame  trap  matrix  a=  However,  in  matrices  of  heterogeneous 
construction  and  even  in  orderly  wire-screen  matrices,  the  free  flow  area 
becomes  a  rather  nebulous  quantity.  In  the  special  case  of  closely  packed 
wire  screens  the  following  approximation  has  been  employed: 

<x=  (1  -  r^DMl  -  n£)  (3-10) 
Further,  if  n^  =  n^  =  n, 

a=  (1  -  nD)2  (3-11) 


Notice  that  the  free  flow  area  obtained  by  using  cx  lies  in  a  plane 
normal  to  the  main  direction  of  flow.  However,  as  in  the  case  of  tube 
banks  it  may  be  desirable  to  employ  the  absolute  minimum  cross-sectional 
area  when  dealing  with  the  pressure  drop  and  heat  transfer  properties  of 
«  matrices j  this  minimum  may  be  in  an  oblique  plane. 


y~h  Equivalent  Diameter 

In  the  case  of  turbulent  flow  inside  tubes,  it  has  been  found  that 
both  pressure  drop  and  heat  transfer  data  obtained  with  passages  of  many 
different  shapes  are  well  correlated  if  the  so-called  equivalent  diameter 
is  employed  as  characteristic  length  in  the  dimensionless  groups.  The 
equivalent  diameter  is 


D 


e 


1;  *  (cross-sectional  area) 
wetted  perimeter 


(3-12) 


In  the  case  of  laminar  flow,  the  use  of  the  equivalent  diameter  has  not 

been  quite  so  satisfactory  and  a  separate  correlation  is  needed  for  each 

shape.  Nevertheless,  much  data  are  often  presented  in  terms  of  D  . 

0 


WADC  TR  55-13 


57 


For  a  wire  screen  matrix,  Eq.  3-12  yields 


«  ha.A  ita 

=  tro ■  s *■ 


(3-13) 


It  follow  that  if  the  screens  are  of  uniform  mesh  and  closely  packed. 


D 

e 


h(l  -  Dn)2 
ir  n 


(3-lli) 


3-5  Matrix  Heat  Capacity 

The  heat  capacity  of  a  matrix  is  equal  to  the  product  of  its  mass 
and  the  specific  heat  of  the  material.  It  follows  that  the  heat  capacity 
per  unit  volume  of  the  matrix  is 

c"  =  f1  '  (3)  f3  =p,s  (3-15) 


3-6  Pressure  Losses 

The  static  pressure  drop  from  the  inlet  duct  to  the  outlet  duct  on 
either  side  of  the  regenerator  nay  be  treated  as  though  the  matrix  were 
a  recuperator,  particularly  if  the  matrix  is  of  the  flame  trap  type.  Thus, 
the  pressure  drop  is  comprised  of  three  parts?  the  entrance  or  contraction 
loss,  the  core  or  frictional  loss,  and  the  exit  or  expansion  loss.  From 
an  equation  used  by  Kays  (Ref.  16)  and  by  Kays  and  London  (Ref.  17),  we 
may  write 


P1~P2  =  T 


ac 


T  i  . 

[(*0+^  Kex)  +  +  a2)(Tjr"_  1)  +  f  qpj  (3-16) 


Kex  are  respectively  the  contraction  and  expansion  co- 


•2g  n 

where  K  '  and 
c 

efficients  based  on  the  dynamic  pressure  inside  the  corej  T  is  the  ab¬ 
solute  temperature,  subscripts  1  and  2  corresponding  to  the  inlet  and 
exit  conditions,  respectively}  f  is  the  friction  factor  of  the  core, 
also  based  on  the  dynamic  pressure  in  the  core}  D  is  the  equivalent 
diameter  of  the  core  passages}  is  the  bulk  mean  temperature  of  the 
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fluid;  [p]  =  lb/ft2;  -  lb/ft3;  and1  [0]  =  lb/sec  ft2. 

In  most  cases,  the  effects  of  the  entrance  and  exit  are  relatively 
small  conpared  with  the  frictional  losses.  Therefore,  it  usually  is  un¬ 
necessary  to  evaluate  K  and  K  with  great  accuracy;  values  in  the 
literature  for  reasonably  similar  situations  may  be  employed.  For 
example,  Kays  (Ref.  16)  has  given  values  for  some  tubular  constructions. 
In  the  case  of  wire  screen  matrices  of  reasonable  length  only  the  last 
two  terms  in  the  brackets  of  Eq.  3-16  need  be  considered  because  the 
other  is,  relatively,  rather  small. 

The  friction  coefficients  f  for  laminar  flow  in  matrices  have  been 
correlated  by  means  of  the  typical  expression. 


where 


G  D 


(3-17) 


(3-18) 


3-6.1  Friction  Factors  for  Flame  Trap  Matrices 

Locke  (Ref.  18)  found  that  the  friction  coefficient  for  smooth 
flame  trap  regenerators  composed  of  round  tubes  lies  within  2  per  cent 
of  the  theoretical  values  obtained  using  C  =  6U. 

.  -V~ 

RomLe  et  al  (Ref.  2k)  performed  pressure  drop  tests  on  two 
flame  trap  matrices,  one  made  of  flat  plates  and  one  of  cylindrical  tubes. 
After  taking  the  entrance  and  exit  losses  into  account,  the  authors  found 
that  the  corresponding  theoretical  values,  namely,  C  =  96  for  the  flat 
plates  and  C  =  6k  for  the  tubes,  were  in  fair  agreement  with  the  experi¬ 
mental  values,  particularly  if  Nj^  «=.  1000. 


The  mass  velocity  G  is  based  on  the  facial  or  superficial  area. 
This  is  to  be  distinguished  from  a  mass  velocity  based  on  the  free  or 
through  area  of  flow.  In  our  symbols  the  mass  velocity  based  on  the  free 
area  is  G/a.  A  bracketed  expression  £...^  is  to  be  read,  "the  units 
of  ...  are". 
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Johnson  (Ref.  15)  performed  single-blow  tests  on  matrix  materials 
of  two  shapes  shown  in  Fig.  3-1 •  The  lengths  were  varied  from  0.073  to 
0.292  ft  by  employing  one,  two,  or  four  elements  in  series.  The  range  of 
the  Reynolds  number  was  2 5  to  200.  His  experimental  results  were  corre¬ 
lated  by  means  of  Eq.  3-17  'with  C  =  17.7  and  n  =  0,78  for  the  matrix 
material  of  0.030-in.  spacing  and  C  =  29.8  and  n  =  0,8?  for  the 
material  of  0.021-in.  spacing.  The  maximum  deviation  of  the  experimental 
points  from  the  mean  correlating  lines  was  about  25  per  cent,  but  on  the 
average  it  was  about  +■  15  per  cent. 

Cox  and  Lamb  (Ref.  ii)  enployed  six  flame  trap  matrices  on  a 
model  rotary  disk  type  regenerator.  The  matrices  were  made  with  triangular 
passages  of  a  type  similar  to  the  ones  used  by  Johnson.  The  ratio  L/D 

© 

took  values  from  31.2  to  82.5,  and  the  Reynolds  number  varied  from 
about  50  to  1*00.  The  authors'  correlation  can  be  represented  by 
Eq.  3-17  if  n  =  1  and 

0  nr* 

c  =  1*1*  +  (1.66  -  ct2  )  (3-19) 

© 

In  their  work,  ol  lay  between  0.7  and  0.83* 

3-6.2  Friction  Factors  for  Wire  Screen  Matrices 

Johnson  (Ref.  15)  performed  tests  on  four  wire  screen  matrices. 
Results  of  the  work  and  a  description  of  the  matrices  are  presented  in 
Fig.  3-2.  It  may  be  noticed  that  the  constructions  of  Matrices  A  and  B 
apparently  differ  only  in  the  total  number  of  layers.  The  results  from 
Matrices  C  and  D,  which  have  the  same  void  fraction,  have  been  correlated 
with  a  single  curve.  The  dotted  line  may  be  represented  by  Eq.  3-17  with 
C  =  25  and  n  =  0.5- 

Jakob,  Kezios,  and  Sogin  (Ref.  ll*)  performed  steady-state  tests 
on  a  wire  screen  matrix  of  0.0075-in.  diameter  wire  and  35  meshj  there 
were  79  layers,  and  the  matrix  length  was  1  inch.  The  authors  found  a 
virtually  uniform  value  for  f,  namely,  1*.8,  in  the  range  of  from 

about  900  to  1200. 
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a  =  0.039 
De*  0.00198  ft 
0.030  spacing. 


a  *  0.800 
De-  0.001324  ft 
0.021  spacing 

FIG. 3-1  DIMENSIONS  OF  FLAME  TRAP 
MATERIAL  OF  REFERENCE  15 
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L 

(ft) 
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(H) 

a 

A 

30 

0.0105 

100 

0.175 

1.66 

0.469 

B 

30 

0.0105 

150 

0.262 

3.57 

0.469 

C 

20 

0.009 

Ipo 

0.150 

1.78 

0.672 

0 

40 

0.0045 

98 

0.094 

1.65 

0.672 

FIG.  3-2  FRICTION  FACTORS  OF  WIRE 
SCREEN  MATRICES  (REF.  15) 
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Ronde  et  al  (Ref.  2k)  observed  that  the  pressure-loss  coefficient, 
in  addition  to  depending  upon  the  Reynolds  number,  depends  upon  the  follow^* 
ing  three  items: 

(1)  The  area  ratio  a.  For  exanple,  the  pressure-loss  coefficient 
was  found  to  double  its  value  when  a  was  reduced  by  half, 

(2)  The  spacing  bertween  screens,  that  is,  the  nunfcer  of  screens  per 
unit  length.  For  exanple,  the  pressure-loss  coefficient  for  an  assembly 
of  unspaced  screens  decreased  1 5  per  cent  when  the  screens  were  separated 
by  a  free  space  of  1/32  inch. 

(3)  The  total  number  of  screens  in  the  assembly.  The  pressure-loss 
coefficient  decreased  or  increased  as  the  nunber  of  screens  increased 
depending,  respectively,  upon  whether  the  screens  were,  closely  or  widely 
spaced. 

More  experimentation  on  configurations  of  wire  screens  are 

needed. 


3-7  Heat  Transfer 

The  coefficient  of  heat  transfer  h,  which  must  be  erployed  in 
evaluating  A  and  (see  Eq.  2-18  and  -19) ,  has  been  found  by 
experiments  based  on  the  single-blow  technique  (e.g.,  Ref.  15,  18,  and  2k) 
and  on  results  from  periodic  regenerators  (e.g..  Ref.  and  20).  Also, 
coefficients  obtained  from  tests  on  certain  recuperator  cores  may  be 
useful  (e.g..  Ref,  17). 

It  has  been  customary  to  correlate  the  coefficients  by  means. of  the 
following  two  typical  equations: 

"Hu  =  °1  >£.  '  (3-20) 

and 

(3-21) 


Nst  =  7^ 
flRe 
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where 


hD 


(3-22) 


Kst  =  TT  (>y) 

P 


and  the  Reynolds  nunber  is  defined  by  Eq.  3-18;  k  is  the  conductivity  of 
the  fluid  evaluated  at  a  mean  tenperature.  Since  the  Stanton  number  is 
the  ratio  of  the  Nusselt  number  to  the  product  of  the  Reynolds  number  and 
the  Prandtl  number  (N^  =  c^yu/k),  the  relationship  between  and 

is 

Cx  =  •  C2  (3-2fc) 

For  air  the  value  of  Nj^  is  0.689  +  0.005  in  the  range  of  tenperature 
from  200  to  800°F. 


In  the  next  two  sections,  data  for  flame  trap  and  wire  screen  matrices 
are  presented.  Regarding  the  relative  merit  of  the  two  types  of  matrices, 
Johnson  remarks  that  the  lowest  weight  for  a  given  performance  is  obtained 
with  a  fine  wire  screen  matrix,  but  the  pressure  drop  is  higher  than  for 
a  smooth  flame  trap  matrix.  The  saving  in  weight  is  more  marked  than  the 
increase  in  pressure  drop;  however,  when  the  partitions  required  in  con¬ 
structing  the  screen  matrices  are  considered,  the  advantage  in  weight  may 
be  reduced. 


3-7*1  Coefficients  of  Heat  Transfer  for  Plane  Trap  Matrices 

Because  the  passages  of  regenerator  matrices  are  small,  the  flow 
is  expected  to  be  laminar.  Glaser  (Ref.  7)  has  shown  that  for  a  circular 
tube  with  linear  tenperature  distribution  (uniform  heating)  in  the  axial 
direction  along  the  wall  and  with  constant  fluid  properties,  the  limiting^- 

^The  limiting  value  of  the  Nusselt  number  is  the  value  toward  which 
the  Nusselt  number  approaches,  by  virtue  of  pure  conduction,  as  the  Rey¬ 
nolds  number  approaches  zero.  Also,  this  is  the  value  that  would  occur 
in  a  very  long  passage  or  in  a  passage  where  the  velocity  and  tenperature 
profile  are  everywhere  fully  developed  (cf.  Eq.  3-25). 


61* 


WADC  TR  55-13 


value  of  the  Nusselt  number  Nv„  „  =  k»36.  Norris  and  Streid  (Ref*  21) 
assured  constant  wall  temperature  and  obtained  ^  =  3,66.  The  ex¬ 

perimental  results  on  cylindrical  passages  by  Locke  (Ref.  18)  fall  between 
the  values  obtained  by  Glaser  and  by  Norris  and  Streid.  The  experimental 
results  of  Romie  et  al  (Ref.  2b)  on  a  circular-tube  matrix  deviate  about 
10  per  cent  from  Glaser’s  theoretical  line  in  the  range  of 
from  7  to  30.  Below  this  range,  down  to  about  3.5,  the  experimental 
data  are  about  20  per  cent  higher  than  the  theoretical  values. 


Clark  and  Kays  (Ref.  3)  performed  calculations  and  steady-state 
experiments  to  study  the  laminar  heat  transfer  in  small  rectangular  pas¬ 
sages  of  various  aspect  ratios.  They  correlated  their  results  by  means 
of  the  typical  equation, 


N„  =  N„  [l  +  C. 
Nu  Nu,  00  |_  ; 


"Re  • 


tt/Dj 


(3-25) 


Both  Njju  ^  and  depend  upon  the  aspect  ratio'*"  Y/X  and  how  the 
tube  wall  is  heated. 


For  a  good  approximation  to  the  regenerator  matrix,  it  may  be 
assumed  that  the  tenperature  distribution  is  practically  linear.  According- 


'3 

N. 


ly,  the  recommended  values  of  N..  and 

Nu,  00 

uniform  heat  input.  The  authors’  values  of 

^  Nu ,  00 

shown  in  Fig.  3-3,  and  their  suggested  values  for  C 
by  the  equation. 


are  those  for  the  case  of 
for  this  case  are 
may  be  represented 


C3  =  0.0025  +  0.020  |  (3-26) 

The  experimental  results  of  Romie  et  al  (Ref.  2k)  on  the  flat  plate  matrix 
(infinite  aspect  ratio)  deviate  about  +  10  per  cent  from  the  values  given 
by  Clark  and  Kays. 

Cox  and  Lamb  (Ref.  Ij)  tested  six  matrices  with  triangular  pas¬ 
sages  in  a  model  rotary  regenerator.  Their  results  for  all  the  matrices 

^See  the  inset  in  the  graph  of  Fig.  3-3. 
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3  LIMITING  NUSSELT  NUMBERS  FOR 
RECTANGULAR  TUBES  WITH  UNIFORM 
HEAT  INPUT 


deviated  by  about  10  per  cent  from  the  correlating  line  represented  by 
Eq.  3-25  with  N^u  ^  =  2.1*  and  =  0.0735.  The  materials  were  like 
those  used  by  Johnson  (Ref.  15).  With  reference  to  Fig.  3-1*  it  may  be 
observed  that  the  passages  were  approximately  isosceles  triangles  with 
included  right  angles.  For  the  equilateral  triangle  with  uniform  heating 
at  the  wall  NNu  ^  =  3.00  (Ref.  3). 

It  has  been  found  inpossible  to  fit  Eq.  3-25  to  all  of  % 

Johnson's  data  for  the  matrices  of  Fig.  3-1,  and  it  is  believed  that  the 

correlation  of  Ref.  4  is  the  more  reliable  one  for  flame  trap  matrices 

with  triangular  passages.  The  author's  original  correlations,  which 

appear  in  the  form  of  Eq.  3-21,  are  tabulated  below.  The  correlations 

represent  the  data  in  the  range  of  from  about  50  to  200.  For 

"Re-  50,  the  Stanton  numbers  appear  to  reach  limiting  values.  Average 

values  of  the  Stanton  number  in  this  range  are  denoted  by  NOJ_  and 

5t,av 

are  presented  in  the  next  to  last  column  of  Table  3-1. 


Table  3-1.  Heat  Transfer  Correlations  for  Flame  Trap  Matrices  (Ref.  15) 


Matrix 

Spacing 

(in.) 

L 

(ft) 

C2 

1-n 

NSt,av 

“NujC 

0.030 

0.0729 

- 

- 

0.01*5 

- 

0.030 

0.11*6 

2.7 

1.0 

0.037 

1.96 

0.030 

0.292 

0.81*5 

0.81 

- 

1.6 

0.021 

0.0729 

- 

- 

0.032 

- 

0.021 

0.292 

0.205 

o.55 

0.021* 

1.5 

Further,  Johnson  remarks  that  at  high  Reynolds  numbers  his  Nus- 
selt  numbers  tend  to  become  constant.  Average  values  of  the  Nusselt  num¬ 
bers  in  this  region  are  denoted  by  c  and  are  presented  in  the  last 
column  of  Table  3-1. 

In  the  so-called  Fraenkl  packing,  two  ribbons  having  corruga¬ 
tions  at  1*5*  to  the  ribbon  edge  are  put  together  so  that  their  corruga¬ 
tions  are  at  right  angles  to  each  other.  Then  the  two  ribbons  are  wound 
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in  a  spiral  to  form  a  disk  of  desired  diameter  and  of  a  height  equal  to 
the  ribbon  width.  One  or  more  such  disks  are  placed  in  series  to  form 
a  matrix  of  desired  length.  Sometimes  spacers  are  used  to  separate  the 
packings.  These  packings  have  been  enployed  mainly  in  regenerators  for 
liquefaction  of  gases.  The  heat  transfer  characteristics  of  Fraenkl  pack¬ 
ings  have  been  investigated  by  Glaser  (Ref.  6)  and  by  Lund  and  Dodge 
(Ref.  20).  Glaser  correlated  his  results  with  the  equation, 

*  jjO.OOTd/Dg)  +  0.65 

"mu  =  ^8.821,(1/0,} +ir.ar  (3_27) 

>7.2  Coefficients  of  Heat  Transfer  for  Wire  Screen  Matrices 

Values  for  h  obtained  by  Johnson  (Ref.  15)  are  in  general 
higher  than  corresponding  values  for  flame  trap  matrices.  Of  course, 
the  pressure  losses  are  correspondingly  higher.  Johnson's  correlations 
can  be  represented  by  means  of  Eq.  3-21.  In  each  case  n  =  0.5.  Values 
of  the  constant  are  presented  in  Table  3-2.  The  range  of  the  Rey¬ 

nolds  number  is  about  25  to  1*00.  Since  some  question  has  been  raised 
regarding  Johnson's  data  for  flame  trap  matrices,  the  same  question  arises 
concerning  his  data  for  wire  screen  matrices.  However,  there  are  no 
other  data  available  on  closely  packed  screens  for  comparison. 


Table  3-2.  Constant  C ^  to  be  Employed  in  Eq.  3-21 
for  Wire  Screen  Matrices  (Ref.  15) 


Matrix 

(cf.  Fig.  3-2) 

Description 

C2 

Mesh 
.  (in.”1) 

Wire  Diameter 

(in.) 

No.  Layers 

A 

30 

0.0105 

100 

0.515 

B 

30 

0.0105 

150 

0.515 

C 

20 

0.009 

100 

1.00 

D 

1*0 

0.001*5 

98 

0.71* 
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Romie  et  al  (Ref.  2U)  employed  two  wire  screen  matrices.  One 
consisted  of  36  layers  of  2h  mesh  screen  made  of  0.00?5-in.  diameter  wire 
and  the  other  of  28  layers  of  16  mesh  screens  made  of  0.013-in.  diameter 
wire.  The  layers  were  spaced  1/32  in.  apart.  Their  Reynolds  number, 
based  on  wire  diameter  as  characteristic  length,  varied  from  about  20 
to  500.  The  authors*  most  reliable  results  are  in  surprisingly  good 
agreement  with  data  for  flow  normal  to  single  wires,  namely, 

“  =  0.32  +  0.1*3  (~±L)  (3-28) 
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Chapter  U:  EFFECT  OF  REGENERATOR  PERFORMANCE 
ON  THE  GAS  TURBINE  POWER  PLANT 

I4-I  Actual  Cycle  of  the  Gas  Turbine  Power  Plant 

The  actual  gas  turbine  cycle  differs  from  the  ideal  cycle  in  several 
nays.  With  reference  to  Fig.  li-1,  which  is  a  temperature-entropy  diagram 
of  an  actual  cycle,  two  importance  differences  may  be  observed: 


il-l.l  Actual  Work  of  Conpressor  and  Turbine 

First,  the  compression  and  expansion  processes  are  isentropic 
in  the  ideal  conpressor  and  turbine  but  not  in  the  actual  ones.  Increases 
of  entropy  occur  in  the  actual  processes.  These  increases  are  represented 
in  the  T,s-diagram  by  the  amounts  of  shifting  of  the  Points  2»  and  I*», 
the  ideal  states,  to  Points  2  and  U,  the  actual  states.  These  changes 
of  entropy  manifest  themselves  as  an  increase  of  the  power  input  to  the 
conpressor  and  a  decrease  of  the  power  output  from  the  turbine.  The  net 
effect  may  be  regarded  as  either  a  decrease  of  the  output  of  a  given 
plant  or  an  increase  of  plant  size  for  a  given  output. 


It  is  custonary  to  express  the  deviations  between  the,  actual 
and  ideal  units  in  terms  of  the  compressor  efficiency, 


hg,  -  hi 


Ui-l) 


and  the  turbine  efficiency. 


-  h3  -  hli 
^t  fau-  hi  . 


( h-2 ) 


I4-I.2  Pressure  Losses 

Second,  there  are  pressure  drops  in  the  low- and  high-pressure 
sides  of  the  cycle  on  account  of  ducting  and  any  heat  exchanger  that  may 
be  present.  Hence,  p^  -=  pg  and  p1<  The  effect  is  to  decrease 

the  power  output  and  the  power  plant  thermal  efficiency. 
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ABSOLUTE  TEMPERATURE  T 


ENTROPY  s 


FIG.  4-1  TEMPERATURE- ENTROPY  PLANE 
FOR  A  REAL  REGENERATIVE 
GAS-TURBINE  CYCLE 
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14-1.3  Leakage 

A  third  difference  between  the  ideal  and  actual  regenerative 
cycle  is  the  leakage  from  the  high-pressure  air  to  the  low-pressure 
exhaust  gas.  Some  is  passed  by  the  seals,  and  some  is  carried  over  by 
the  rotation.  On  account  of  this  leakage,  the  amount  of  fluid  passing 
through,  the  conpressor  is  greater  than  that  passing  through  the  turbine. 
Hence,  a  larger  and  more  powerful  conpressor  is  required  for  a  given  tur¬ 
bine  output.  The  net  result  is  a  decrease  in  the  thermal  efficiency  of 
the  power  plant. 

h-2  Thermal  Efficiency  of  the  Actual  Cycle 

The  power  output  of  the  plant  and  the  thermal  efficiency  of  an  actual 
gas  turbine  cycle  operating  with  practical  conpression  and  expansion  units, 
with  pressure  drop  in  the  intermediate  ducting  and  regenerator,  and  with 
leakage  from  the  air  side  to  the  gas  side  are  developed  in  Appendix  B. 

The  calculations  are  based  on  the  assumptions  that  the  fluid  is  thermally 
and  calorically  perfect  and  that  the  fuel  has  negligible  influence  on  the 
ther modynami c  p rop er tie  s . 


U— 3  Evaluation  of  the  Effects  of  Leakage  and  Pressure  Drop 

Harper  and  Rohsenov/  (Ref.  $)  employed  the  results  of  Appendix  B  or 
their  equivalents,  to  calculate  and  demonstrate  the  effects  of  leakage 
and  pressure  drop  for  the  following  operating  conditions: 


Air  inlet  temperature, 
Maximum  temperature, 
Conpressor  efficiency. 
Turbine  efficiency, 


60°F,  or 
1500 *F,  or 


520 'R 
1960*R 
0.85 
0.88 


The  authors  expressed  their  assumed  values- of  leakage  and  pressure  loss, 
respectively,  in  terms  of  the  leakage  fraction  Auy^  and  the  sum  of 
the  pressure  drop  ratios  ]>]  (Ap/p).  The  quantity  m  is  the  mass  rate 
of  flow  through  the  conpressor;  and  with  reference  to  Fig,  ij-1. 


P2  ~  p3  PU  ~  P1 
P2  Pi 


{hr  3) 
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The  ratio  of  the  first  term  of  Eq.  U—3  to  the  second,  which  would  be  de¬ 
noted  by  a^/cc^  in  Appendix  B,  was  not  given  by  Harper  and  Rohsenow  but 
appears  to  have  been  taken  as  3/7  in  one  part  of  their  calculations* 

The  authors  calculated  the  maximum  thermal  efficiency  and 

the  corresponding  (optimum)  pressure  ratios  for  the  cases  that  An/m  =  0, 
0.05,  and  0.10  and  that  (Ap/p)  =  0.05,  0.10,  and  0.15.  Results 
of  their  study  are  discussed  in  the  next  two  sections. 


lj-3.1  Effect  of  Leakage 

As  has  already  been  remarked  in  Section  2-13,  the  influence  of 
leakage  on  the  effectiveness  of  the  heat  exchanger  is  small.  The  influence 
of  leakage  on  the  poirer  plant  performance  is  shown  in  Fig.  U-2a  and  -2b. 
With  regard  to  increasing  values  of  An/™*  the  following  changes  rray  be 
observed : 

(1)  The  maximum  attainable  thermal  efficiency  decreases.  This  is 
mainly  due  to  the  decrease  in  turbine  output,  for  the  mass  flow  in  the 
turbine  is  less’  than  the  mass  flow  at  the  intake  of  the  power  plant. 

(2)  The  optimum  pressure  ratio  r  ,  decreases.  The  effect  is 

P,°Pt 

most  pronounced  at  lowest  values  of  the  regenerator  effectiveness. 

(3)  The  work  output  per  pound  of  air  entering  the  power  plant  de¬ 
creases.  This  is  due  to  the  decrease  of  both  the  mass  flow  in  the  tur¬ 
bine  and  the  optimum  compression  ratio. 


h-3 . 2  Effect  of  Pressure  Drop 

The  effect  of  pressure  drop  is  shown  in  Fig.  U-3a  and  -3b.  With 
regard  to  increasing  values  of  2  (Ap/p),  the  following  changes  may  be 
observed: 

(1)  The  maximum  attainable  thermal  efficiency  decreases.  This  is 
due  to  the  increase  in  the  turbine  ejqjansion  ratio. 

(2)  The  optimum  pressure  ratio  increases. 

(3)  The  work  output  per  pound  of  air  decreases. 
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TURBINE  POWER  PLANT  FOR  PRESSURE  DROP  AP,  *0.10 


l*-li  Effects  of  Regenerator  Speedy  Length,  and  Mass  Velocity 

Harper  and  Rohsenow  (Ref.  9),  employing  the  results  of  Saunders  and 
Smoleniec  (Ref.  25),  also  showed  how  the  rotational  speed  N,  regenerator 
length  L,  and  the  mass  velocity  m"  affect  the  power  plant.  In  con¬ 
sidering  any  one  of  these  three  variables,  the  other  two  were  regarded 
as  constant.  Also,  the  net  output  of  the  gas  turbine  power  plant  was 
maintained  constant  under  the  conditions  specified  in  Section  ij-3.  A 
balanced  regenerator  with  A  =  5  and  ft  =  k  was  taken  as  a  basis 
for  reference.  The  results  of  the  authors'  discussion  are  summarized 
in  the  next  three  sections. 

li— u.l  Effect  of  Regenerator  Speed 

An  increase  of  regenerator  speed  decreases  I  proportionately 
but  leaves  A  unaltered.  The  increase  of  speed  increases  the  carry-over 
losses.  It  also  influences  the  leakage  through  the  seals  insofar  as  the 
regenerator  speed  affects  the  optimum  pressure  ratio.  Regarding  the 
pressure  losses,  Harper  and  Rohsenow  took  ^  (Ap/p)  to  be  a  fixed  value. 

The  results  of  their  calculations  are  shown  in  Fig.  k-k»  It 
may  be  observed  that  the  maximum  cycle  efficiency  reaches  a  maximum  at  a 
certain  speed;  then  it  decreases  on  account  of  the  increasing  carry-over 
losses.  The  matrix  relative  volume  V  increases  since  the  work  per 
pound  of  air  W  decreases  while  total  power  output  remains  constant. 

The  optimum  pressure  ratio  r^  ^  decreases  as  the  speed  increases. 

li-U.2  Effect  of  Regenerator  Length 

In  this  case  the  reduced  length  A  changes  while  the  reduced 
period  1  remains  constant.  The  results  of  calculations  for  a  fixed 
value  of  An/®  are  shown  in  Fig.  lj-5. 

The  pressure  drop  (Ap/p)  increases  with  the  matrix  length. 
Again,  the  naximum  thernal  efficiency  reaches  an  optimum  value;  then  it 
decreases  on  account  of  the  rising  pressure  drop.  Since  the  work  per 
pound  of  air  entering  the  turbine  (denoted  by  W)  decreases  as  the  matrix 
length  increases,  the  matrix  facial  area  as  well  as  length  increase  to 
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FIG. 4-4  EFFECT  OF  ROTOR  SPEED  ON 
GAS  TURBINE  POWER  PLANT 
PERFORMANCE  AT  DESIGN  - 
POINT  CONDITIONS 
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FIG.  4-5  EFFECT  OF  MATRIX  LENGTH  ON 
GAS  TURBINE  POWER  PLANT 
PERFORMANCE  AT  DESIGN  -  POINT 
CONDITIONS 
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.  •  I 

-  .  I 

) 

provide  a  given  fixed  total'  power  output.  For  this  reason  the  relative  ; 

voluaB  V  increases  more  rapidly  than  just  linearly. 

h-h.3  Effect  of  Mass  Velocity 

A  change  of  the  rate  of  flow  nay  affect  both  the  reduced  length 
and  the  reduced  period*  because  it  nay  be  accompanied  by  a  change  in  the 
value  of  the  coefficient  of  heat  transfer  h.  For  example*  if  it  is 
assumed  that  the  flow  through  the  matrix  is  such  that 

h  ~  (m")^2  (Wi) 


0H>) 

and 

It2  ~  a"  (h-6) 

Assuming  that  the  pressure  drop  increases  linearly  with  m"  and  that 
A  Wm  remains  constant*  Harper  and  Rohsenow  presented  the  results  shown 
in  Fig.  h-6.  It  can  be  seen  that  W  reaches  a  maximum  value  after  which 
it  decreases  on  account  of  the  rising  pressure  loss. 


then 

A2  ~  (m")”1 


\ 
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(PROPORTIONAL  TO  THE  RATE  OF  MASS  FLOW  OF  FUJD) 

FIG.  4-6  EFFECT  OF  RATE  OF  MASS  FLOW 
ON  GAS  TURBINE  POWER  PLANT 
PERFORMANCE  AT  DESIGN— POINT 
CONDITIONS 
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Chapter  5:  A  DESIGN  PROCEDURE 


5-1  Preliminary  Remarks 

As  mentioned  in  the  introductory  chapter,  this  manual  deals  ire  inly 
•with  the  thermal  design.  However,  many  features  of  the  mechanical  design 
dictate  the  thermal  performance;  and  limitations  of  weight,  space,  and 
location,  in  contributing  to  the  over-all  problem,  may  very  well  deter¬ 
mine  whether  the  heat  exchanger  is  a  practicable  matter  in  a  given  appli¬ 
cation. 

The  designer  must  use  his  science  and  art  to  build  a  regenerator  that 
lies  within  the  allowable  mechanical  limitations  and  that  provides  a 
satisfactory  over-all  performance  of  the  power  plant.  He  faces  this  task 
with  a  system  or  group  of  quantities  whose  values  are  given.  Some  of 
these  values  have  had  to  be  determined  from  the  performance  requirements, 
some  by  auxiliary  calculations,  and  others  by  judgment,  which  improves 
with  experience.  Generally,  his  design  procedure  depends  upon  the  par¬ 
ticular  group  of  values  he  has  to  start  with. 

In  the  following  sections  a  rather  special  problem  is  treated.  A 
statement  of  the  problem  appears  in  the  next  section.  Then  the  analysis 
is  briefly  discussed.  This  is  followed  by  a  solution  with  some  sample 
calculations.  The  solution  illustrates  the  application  of  much  of  the 
material  in  Chapters  2  and  3*  The  results  are  not  intended,  however,  to 
be  definitive  without  further  development  and  experimentation,  particularly 
with  regard  to  the  pressure  drop  and  heat  transfer  characteristics  of  the 
type  of  matrix  selected.  It  is  hoped  that  the  procedure  presented  here 
may  suggest  how  other  groups  of  starting  quantities  may  be  handled  and, 
consequently,  how  other  design  procedures  may  be  devised. 

5-2  Statement  of  the  Problem 

life  suppose  that  the  designer  has  been  given  the  quantities  in  Table  5-1. 
As  can  be  seen,  his  regenerator  is  to  achieve  an  effectiveness  of  90  per 
cent,  suppose  that  a  drum  type  regenerator  is  to  be  designed.  For  any 
matrix  construction  that  he  selects  the  designer  must  calculate  the  quanti¬ 
ties  listed  below* 
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L  the  matrix  length,  that  is,  the  thickness  of  the  drum;  ft 
the  mean  or  pitch  diameter  of  the  drum;  ft 
\  the  axial  length  of  the  drum;  ft 
2  it  K  the  angle  subtended  by  the  air  side 

ft 

2  itK  the  angle  subtended  by  the  gas  side 

O 

N  the  speed  of  rotation;  revolutions  per  second 


Table  $-1.  Given  Design  and  Performance  Conditions 
for  the  Illustrative  Problem 


Symbol* 

Quantity 

Value 

Units 

P 

power  at  the  shaft 

2500 

hp 

1c 

conpressor  efficiency 

0.85 

— 

*1  confc 

confcustion  efficiency 

0.95 

— 

'It 

turbine  efficiency 

0.85 

— 

^Reg 

required  regenerator  effectiveness 

0.90 

— 

Ti 

tenperature  at  conpressor  intake 

520 

*R 

T2 

tenperature  at  regenerator  inlet  (air  side) 

815 

•r 

T3 

tenperature  at  throttle 

1810 

*R 

V 

tenperature  at  regenerator  inlet  (gas  side) 

1376 

•r 

P2 

pressure  at  regenerator  inlet  (air  side) 

58.8 

lb/sq  in. 

Pi* 

pressure  at  regenerator  inlet  (gas  side) 

15.1*5 

lb/sq  in. 

Apa 

allowable  pressure  drop  on  air  side  . 

0.88 

lb/sq  in. 

*Pg 

allowable  pressure  drop  on  gas  side 

0.75 

lb/sq  in. 

m 

a 

rate  of  flow  on  air  side 

11*0,000 

lb/hr 

m 

g 

rate  of  flow  on  gas  side 

11*1,000 

lb/hr 

Subscripts  1,  ...,  1*  refer  to  the  statepoints  in  the  cycle  diagram 
of  Fig.  k-1.  All  tenperatures  and  pressures  are  absolute. 
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5-3  Analysis  of  the  Solution 

After  a  satrix  is  selected  the  basic  geometric  and  thermal  properties 
can  he  evaluated  in  accordance  with  Sections  3-1  to  -5.  General  relation¬ 
ships  regarding  pressure  drop  and  heat  transfer  data  can  be  selected  from 
Sections  3-6  and  -7. 


The  pressure  drop  depends  upon  the  matrix  length  L  and  the  mass 

velocities^-  G  and  G  .  This  gives  two  relations,  one  between  L  and 
a  g  6 

G_  and  the  other  between  L  and  G  . 

ft  g 


The  heat  transfer  coefficients  h  and  h  depend,  respectively, 

& 

Therefore,  the  reduced  lengths  A  and  A^  can 


upon  G  and  G 


g 


g 


be  expressed  in  terms  of  G&,  G  ,  and  L.  Enploying  the  first  two  re¬ 


lations,  we  may  eliminate  G  and 


_  G  to  express  and  A  in 

ft  g  ft  g 

terms  of  L.  Consequently,  the  harmonic  mean  reduced  length  A  can  be 

expressed  in  terms  of  L. 


m 


The  allowable  carry-over  losses  Am  from  the  cold-air  side  depends 
upon  the  speed  of  rotation  and  the  volume  of  the  matrix.  It  will  be  seen 
that  a  known  value  of  Am  /m  provides  a  relation  among  the  quantities 
N,  G  ,  and  L.  Also,  the  utilization  factor  depends  upon  N,  G  >  and  L. 
Observing  that  the  fractions  Kr  and  are  sinply  related  to  Gft  and 

G  ,  we  can  calculate  a  unique  value  of  the  utilization  factor  U.  Hence, 

o 

Affl  can  be  found  from  Fig.  2-8,  and  all  required  quantities  can  be  cal¬ 
culated. 


5-U  Selection  of  the  Matrix 

The  uatrix  of  our  regenerator  is  to  be  made  of  wire  screen.  It  has 
been  suggested  that  the  material  forming  the  matrix  is  limited  by  burnout 
due  to  fuel  left  in  the  exhaust  gases  and  that  the  lower  limit  in  the  size 

p 

of  wires  (or  plates)  may  be  in  the  order  of  0.008  inch. 

^It  has  been  found  convenient  to  let  fm"!  =  lb/hr  ft^  and  FgI  = 
lb/ sec  ft^j  hence,  G=  3600  m".  The  brackets  here>  as  elsewhere, 

are  to  be  read,  "the  units  of  ...  are". 

^See  the  reference  to  the  work  of  Bowden  and  Hryniszak  at  the  end  of 
the  Bibliography, 
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i 


Johnson  (Ref.  1 5)  found  that  copper  gauze,  glass  wool,  and  steel  wool 
are  unsuitable  for  high  tenperatures.  Thin  strips  of  a  70:30  copper-nickel 
alloy  and  wire  gauzes  of  brass,  of  mild  steel,  and  of  stainless  steel 
showed  no  deterioration  or  corrosion  after  tests  at  about  900*F. 

For  purposes  of  the  illustrative  calculations,  the  gauze  is  assumed 
to  be  made  of  stainless-steel  wire  0,009  in.  in  diameter.  It  is  supposed 
that  the  screens  are  of  20  mesh  and  are  firmly  packed.  The  basic  geometric 
and  thermal  properties  have  been  calculated  by  neans  of  the  relationships 
in  Chapter  3.  The  results  are  presented  in  Table  5-2. 


Table  5-2.  Basic  Geometric  and  Thermal  Properties 
of  the  Assumed  Matrix  Material 


Symbol 

Quantity 

Value 

Units 

D 

wire  diameter 

0.009 

0.00075 

in. 

ft 

n 

mesh 

20 

in.”1 

P 

void  fraction 

0.86 

- 

1  -  (3 

solid  fraction 

o.il* 

- 

SI« 

heat  transfer  area  per  unit  volume 

755 

ft”1 

a. 

ratio  of  free-flow  area  to  facial  area 

0.67 

- 

?s 

density  of  the  solid 

U87 

lb/ft3 

°p,s 

specific  heat  of  the  solid 

0.138 

B/lb  F 

Cm 

heat  capacity  of  the  matrix  per  unit  volume 

9k 

B/f t3  F 

D 

e 

equivalent  diameter 

0.00355 

ft 

5-5  Calculation  of  the  Mean  Inlet  and  Outlet  Tenperatures  on  the  Cold 
and  Hot  Sides  of  the  Regenerator 

By  definition  of  the  regenerator-  effectiveness,  the  mean  outlet  tem¬ 
perature  on  the  cold  side  of  the  regenerator  is 

TA  =  T2  +  WT1,  "  V  =  1320#R 
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Similarly,  the  mean  outlet  tenperature  on  the  hot  side  is,  practically, 

TB  =  -  W*k  -  V  =  871‘a 

Table  5-3  is  a  resume  of  the  fluid  temperatures  at  the  regenerator 
inlets  and  outlets.  Symbols  employed  in  the  thermodynamic  cycle  of  Fig.  b-1 
are  related  to  the  corresponding  symbols  employed  in  Chapter  2.  Thus, 

T2  =  8l5*R  while  ©&  1  =  Tg  -  h60  =  355*F,  and  so  forth. 


Table  5-3 •  Resume  of  Temperatures  at  the  Regenerator  Inlets  and  Outlets 


Cold  Air 

Hot  Gas 

in 

out 

in 

out 

T2  =  8l5*R 

,  =  355* F 

a,l 

T.  =  1320*R 

A 

0  .  =  860*F 
a,2 

=  1376*R 

©  ,  =  916* F 
g>l 

Tfi  =  871*R 
9  =  l4ll*F 

5-6  Evaluation  of  the  Mean  Tenperature  of  Each  Fluid 

In  accordance  -with  a  recommendation  of  Cox  and  Lamb  (Ref.  4),  the 
fluid  properties  with  the  exception  of  the  specific  weight,  are  evaluated 
at  a  tenperature  which  is  0.6  of  the  way  between  the  average  wall  tem¬ 
perature  and  the  bulk  mean  fluid  tenperature.  This  rule  is  used  here  even 
though  Cox  and  Lamb  worked  only  on  flame  trap  matrices  and  Johnson  corre¬ 
lated  his  data  with  a  somewhat  different  definition  of  the  mean  tenperature. 


The 


bulk  mean  air  temperature  is  0  =  (©  _  +  ©  _)/2  =  608*Fj 

&  jD  &  jX  B  y  d. 


the  bulk  gas  tenperature  is  ©^  ^  =  (©^  ^  +®g  ?)  f  ^  ~  ^6it*F. 

The  average  wall  tenperature  t  is  assumed  to  be  midway  between 

©  ,  and  ©  ,  .  Accordingly,  in  the  illustrative  example,  t  =  636*F. 
gjD  s,m 

Therefore,  the.  tenperature  to  be  used  for  evaluating  the  fluid  proper¬ 
ties  (y*.,  k,  and  c^)  on  the  air  side  is 

©  =  t  -  0.6(t  -  ©  J  =  625*F 

a,m  s,m  '  s,m  a,b' 
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•while  on  the  gas  side  it  is 


©  =  t  +  0.6(0  .  -  t  )  =  6h7*F 

g,m  s,m  g,b  s,m 


The  specific  weights  are  evaluated  at  the  bulk  ms an  stream  tenpera- 
tures  under  the  assumptions  that  both  the  cold  and  the  hot  fluids  are 
air  and  that  air  is  a  thermally  perfect  gas  (gas  constant  R  =  53.3  ft/F). 
The  mean  pressures  on  the  air  and  gas  sides  are  8I1IO  and  2170  lb/ft^, 
respectively.  Table  5-ii  is  a  resume  of  the  fluid  property  values  that  are 
employed  in  the  calculations  which  follow. 


Table  5-lu  Values  of  the  Fluid  Properties  in  the  Illustrative  Example 


Quantity 

Symbol 

fluid 

Units 

Cold  Air  , 

Hot  Gas 

T 

0.168 

0.0362 

lb/ft3 

c 

P 

0.251 

0.252 

B/lb  F 

10^  •  ^ 

20.  k 

20.6 

lb/sec  ft 

k 

0.0270 

0.0271 

B/hr  ft  F 

5-7  Relationships  Derived  from  the  Correlation  of  Pressure  Drop  Data 
From  Section  3-6.2,  the  pressure  drop  of  the  natrix  is,  in  a  first 
approximation, 


=  t  ■  IT 

e 


where,  in  the  same  approximation. 


(5-1) 


r  =  — — 


(5-2) 


e 

r 


Eliminating  f  and  introducing  the  numerical  values  of  g,  a,  and  D  , 

© 
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(5-3) 


0.5 

Ap  =  3200  -£■ —  •  L  •  ff 

r 

where  [Ap]  =  lb/ft2,  [^]  =  lb/sec  ft,  [y]  =  lb/ft3,  [L]  =  ft,  and 
[G]  =  lb/sec  ft2.  With  reference  to  Table  5-1,  the  allowable  pressure 
drops  on  the  air  and  gas  sides  are  127  and  108  lb/ft2,  respectively. 
Therefore,  Eq.  5-3  yields  the  system  of  equations, 


G  =  1.1*75  L“2/3 

a 

(5-10 

Gg  S  0.1a8  L"2/3 

(5-5) 

5-8  Calculation  of  the  Angle  Subtended  by  the  Ducts 

Since  the  leakage  fraction  and  the  ratio  of  fuel  to  air  are  small. 

9  l 

it  is  assumed  that  U&  =  Ug  =  U,  in  accordance  with  Eq.  2-25.  By  defi¬ 
nition. 


KGc 


U  = 


Cm  LN 


Since  c  «  c  ,  it  follows  that 
P>a  P»g 


(5-6) 


K  G 


a  a 


K  G 
g  g 


(5-7) 


Further, 


K  +  K  =  1 
a  g 


K 


sis 


(5-8) 


where  Kgls  is  the  fraction  of  the  circle  allowed  for  seals,  duct  walls, 
etc.  Assuming  that  Kglg  =  0.06  (corresponding  to  an  angle  of  about  22*) 
and  eliminating  Gft,  G  ,  and  L  from  Eq.  5-ii,  -5,  and  -7,  we  obtain  the 
set  of  simultaneous  equations 


Ka  +  Kg  =  0.9U 

K  -  0.285  K  =  0 
a  g 


(5-9) 

(5-10) 
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The  solution  is  K  =  0.21  and  K  =  0.73,  corresponding  to  angles  of 

a  g 

about  75*  and  263*,  respectively.  Adjustments  may  have  to  be  made  at 
a  later  stage  of  the  design. 


5-9  Relations  Derived  from  the  Correlation  of  Heat  Transfer  Data 

From  Section  3-7.2,  the  correlation  of  the  heat  transfer  data  for  the 
selected  matrix  is 


cc 


1.00 


3600  G  c 


nn 

\<x  • 


(5-11) 


r 


Substituting  the  numerical  values  given  in  the  preceding  sections  for  ccj 


r- 


and  0^,  we  obtain  the  system  of  equations, 


h  =  83.5  G 

a  a 


(5-12) 


h  =  8 It. 3  G 
g  g 


where  =  B/hr  ft2  F  and  QjQ  =  lb/ sec  ft2, 
-with  the  expressions  of  Eq.  5-U  and  -5, 


h  =  101  L 

a 


-1/3 


h  =  5U.6  L 

o 


-1/3 


(5-13) 


Eliminating  G  and  G 
a  g 


(5-lW 

(5-15) 


By  definition,  the  reduced  length  is 


A  = 


h  SmL 
3600  G  c 


(5-16) 


Since  h  and  G  are  known  in  terms  of  L,  the  reduced  lengths  may  also 
be  expressed  in  terms  of  L.  Thus,  enploying  Eq.  5— It,  -5,  -7,  -8,  and  -11, 


Aa  =  57.lt  L 


V3 


A  =  108  L 

& 


V3 


(5-17) 

(5-18) 
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where  [[L^  =  ft.  Further,  the  harmonic  mean  reduced  length  (cf.  Section 
2-9.1)  is 


Am  =  75.3  L1^3  (5-19) 

The  regenerator  could  now  be  uniquely  determined  if  another  relation 
between  A  and  L  were  available.  In  the  next  section  an  indirect 

SI 

relationship  is  obtained:  the  utilization  factor  U  is  evaluated  in 
terms  of  the  allowable  carry-over  loss.  Its  value  is  all  that  is  needed, 
because  Affi,  U,  and  A  Fleg  are  re^ate^  by  the  theory  of  the  balanced  re¬ 
generator. 


5-10  Determination  of  the  Utilization  Factor  in  Terms  of  the  Carry-Over 
During  each  revolution  of  the  matrix,  cold  air  is  carried  over  to 
the  gas  side  and  hot  gas  to  the  air  side.  For  purposes  of  illustration 
it  will  be  assumed  that  only  the  carry-over  of  the  cold  air  need  be  con¬ 
sidered,  because  its  density  is  much  greater  than  that  of  the  hot  gas. 
Also,  it  will  be  assumed  that  the  leakage  at  the  seals  is  so  small  that 
only  the  carry-over  is  significant;  in  practice  this  may  not  be  the  case 

and  an  estimate  of  the  leakage  at  the  seals  would  have  to  be  made  during 

1 

some  stage  of  the  design  or  development. 

The  rate  of  loss  of  cold  air  to  the  hot  side  by  virtue  of  the  carry¬ 
over  equals  the  speed  of  rotation  times  the  weight  of  the  air  contained 
in  the  total  volume  of  the  void  space  in  the  matrix.  Thus, 

=  3600  »  .  Vtot  •  p  •  u,  (5-20) 

where  V,  .  is  the  total  volume  of  the  matrix.  In  terms  of  A  ,  to  a 
tot  a' 

good  approxiaation, 

7t0t  =  =  V(1  +  (5-a) 

a 

■^An  analysis  of  the  seals  may  be  found  in  the  work  of  Bowden  and 
Hryniszak.  See  the  Additional  References  at  the  end  of  the  Bibliography. 
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This  allows  for  partitioning  the  matrix.  Substituting  this  expression 

for  V.  .  into  Eq.  5-20,  observing  that  m  =  36OO  G  A  ,  and  rearranging 
"tot  ft  ft  ft 

the  resultant  equation. 


G 

a 

LN 


K 


*  Vs'1**6) 


Am 

(-—) 


m_ 


(5-22) 


The  influence  of  Am  /m  on  the  power  plant  efficiency  could  be  checked 
by  enploying  the  material  in  Chapter  i*.  This  is  not  done  here,  and  for 
the  sake  of  brevity  in  the  illustrative  exanple  it  is  assumed  that  a  satis¬ 
factory  value  of  Am  /m  has  been  found  to  be  1  per  cent.  From  Fig.  U-2 

ft  ft 

it  can  be  seen  that  if  the  leakage  through  the  seals  is  also  quite  small, 
no  correction  is  necessary  for  the  regenerator  effectiveness. 


It  follows  that 

rnr  =  6lx 

where  =  lb/sec  ft  ,  [V]  =  ft,  and  DC  =  sec”  .  The  right-hand 

memb.er  of  Eq.  5-6  is  now  evaluated  with  the  result  that  U  =  0.036. 


5-11  Results  of  the  Sanple  Calculations 

Entering  the  theoretical  curves  of  Fig.  2-8  for  the  balanced  regenera¬ 
tor  with  =  °«90  and  U  =  0.036,  we  find^  Am  =  18.  (It  may  be  ob¬ 

served  that  this  is  exactly  the  value  that  would  have  been  obtained  by 
employing  Eq.  2-72,  the  reason  being  that  the  utilization  factor  in  this 
illustrative  example  is  very  small.)  Consequently,  from  Eq.  5-15, 

L  =  O.3U2  ft  or  about  U.l  inches.  Having  found  the  value  of  L,  we 
may  readily  obtain  all  other  quantities  of  interest.  The  values  of  many 
of  them  are  listed  in  Table  5-5. 
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Table  5-5.  Results  of  the  Calculations 


Symbol 

Air  Side 

Gas  Side 

Units 

U 

0.036 

0.036 

— 

L  ' 

0.3U2  ■ 

0.31*2 

ft 

G 

3.02 

0.855 

lb/sec  ft2 

h 

11*5 

?8 

B/hr  ft2  F 

N&>  =  ®e/(A,C° 

780 

222 

— 

Nst  =  h-a/(G  cp) 

0.0356 

0.0658 

— 

f 

0.86 

1.61 

— 

A 

13.7 

25.8 

— 

U 

0.1*91* 

0.929 

— 

K 

0.21 

0.73 

— 

N 

0.137 

0.137 

rps 

N 

8.22 

8.22 

rpm 

Prom  the  relationship  between  the  rate  of  flow  m  and  the  mass 

2 

velocity  G,  it  is  found  that  A  and  A  are  12.9  and  1*5.8  ft  ,  re- 

®  £ 

spectively.  The  area  allowed  for  the  seals  being  about  6  per  cent  of  the 

2 

total  area,  the  pitch  surface  of  the  drum  has  a  total  area  of  62.5  ft 
and  the  volume  of  the  drum  is  21.1*  ft^.  >  Assuming  that  the  length  •  * 

of  the  drum  is  1.5  times  its  pitch  diameter  D^,  one  finds  that 
Dj  =  3.65  ft  and  \  =  5.1*6  ft.  The  matrix  screens  would  weigh  about 
one  ton. 

5-12  Final  Remarks 

It  has  been  noted  that  an  assumption  underlying  the  calculation  pro¬ 
cedure  is  that  the  pressure  drop  and  the  heat  transfer  correlations  of 
the  selected  matrix  material  can  be  extended  beyond  the  range  of  Johnson’s 
tests.  This  would  require  an  experimental  check.  In  particular,  the 
assumed  friction  factor  on  the  cold  side  appears  to  be  considerably  lower 
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than  a  true  extrapolation  of  the  data  indicates.  The  sanple  calculation 
could  be  inproved  by  taking  this  difference  into  account* 

The  value  0.036  for  the  utilization  factor  is  rather  low.  This 
indicates  that  the  carry-over  losses  could  be  reduced  by  lowering  the 
matrix  speed  without  seriously  inpairing  the  effectiveness  of  the  cal¬ 
culated  matrix  or  without  having  to  increase  the  matrix  length  consider¬ 
ably  to  retain  the  same  effectiveness.  In  general,  it  is  expected  that  a 

t 

regenerator  could  be  found  to  operate  satisfactorily  with  a  utilization 
factor  in  the  order  of  magnitude  0.3. 

The  influence  of  conduction  along  the  x-direction  of  the  matrix  can 

be  calculated  by  means  of  the  chart  in  Fig.  2-9.  In  the  illustrative 

exanple,  the  conductivity  of  the  stainless  steel  is  k  =  15  B/hr  ft  F, 

s 

and  so  »  0.008.  Therefore,  the  influence  of  the  conduction 

max 

in  the  x-direction  is  to  decrease  the  effectiveness  less  than  1  per  cent. 

The  designer  may  have  to  try  a  few  matrix  materials  before  he  is 
satisfied  that  he  has  the  regenerator  of  least  weight  and  volume  for 
his  application.  Having  found  a  matrix  to  satisfy  the  thermal  require¬ 
ments,  the  designer  may  proceed  to  the  details  of  the  mechanical  design. 
Finally,  he  can  return  to  his  calculations  of  the  thermal  performance, 
inproving  his  estimates  and  predicting  the  performances  at  other  speeds 
and  other  loads. 
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INFLUENCE  OF  CURVATURE  IN  DRUM  TYPE  MATRICES 
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Appendix  At  INFLUENCE  OF  CURVATURE  IN  DRUM  TYPE  MATRICES 


A-l  General  Remarks 

In  this  appendix  it  is  shown  that  the  perforirance  of  the  drum  type 
regenerator  can  be  obtained  directly  from  the  theory  of  the  flat  regenera- 
,  tor  (cf.  Chapter  2).  The  demonstration  is  based  on  the  assumptions  of 
Section  2-1  and  is  carried  out  using  the  line  of  reasoning  in  Section  2-2. 

A-2  Derivation  of  the  Differential  Equations 

A  drum  type  matrix  of  unit  length  along  the  axis  is  considered.  The 
angular  speed  is  to  .  Heat  balances  are  made  on  an  element  of  space  shown 
in  Fig.  A-l.  Since  the  mass  velocity  mn  varies  with  the  radius,  it 
■will  be  convenient  to  employ  the  axial  interstitial  velocity  u,  because 
u  •  r  is  a  constant  by  continuity.  Also,  the  matrix  velocity  v  is  a 
function  of  the  radial  distance,  namely,  v  =  ud  •  r. 

It  is  sufficient  to  consider  only  one  of  the  fluids;,  and.  it  may  be 
supposed  that  its  channel  subtends  an  angle  K  •  2iv  radians.  The  solid 
teEperature  is  t  =  t(r,0)  and  the  fluid  temperature  0  =  O(r,0) . 

A-2 .1  Thermal  Energy  Conveyed  by  the  Fluid 

The  energy  transported  per  unit  time  by  the  fluid  into  the  ele¬ 
ment  of  space  is 

qj.  =  U  •  (r  60  .  1)  p  cp  •  0  (A-l) 

A-2. 2  Thermal  Energy  Conveyed  by  the  Solid 

Similarly,  the  thermal  energy  transported  per  unit  time  into  the 
element  of  space  by  the  solid  material  by  virtue  of  its  enthalpy  is 

^  l  =  u>  r  •  (5r  •  1)  •  Cw  •  t  (A-2) 
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A-2.3  Thermal  Energy  in  the  Carry-Over 

The  thermal  energy  of  the  fluid  carried  over  per  unit  time  is 

q02  =  ^u5r*(p«5r*l)  •  cp  •  ©  (A-3) 


A-2.li  Balance  of  the  Transported  Thermal  Energy 

Applying  the  First  Law  of  Thermodynamics  to  the  elemental  volume 
of  space  and  observing  that  no  heat  is  stored, 


5r  .  ^  ®  -  0 

dr  50 


(A— It) 


Substituting  from  Eq.  A-l,  -2,  and  -3  into  Eq.  A-ii  and  observing  that 
p,  c  ,  and  u  •  r  are  constants, 


d© 


u  p(3  c  — 
V  '  P  dr 


+  C'"  oo 


LP  c 


d©, 

P  dj? 


=  0 


(A-5) 


A-2.5  Heat  Exchange  in  the  Elemental  Volume  of  Space 

The  rate  of  change  of  the  enthalpy  of  the  matrix  material  is 
equal  to  the  rate  of  heat  transferred  by  convection  in  the  same  space 
element.  Accordingly, 


d 


(Cm  ip  r 


t  •  5r)  50  =  hS'" 


(©  -  t)  •  r  6r  60 


(A-6) 


or 


at  _  hsw  . 

50  -  "C»  "  t) 


(A— 7) 


Eliminating  dt/90  from  Eq.  A-5  and  -7, 


uePcp^  +  fPwop^=hS"(t-e) 


(A-6) 


It  is  convenient  to  eirploy  the  rate  of  fluid  flow  per  unit  axial  length. 


m 


'  =  u  •  (K  •  2 v  r  •  p) 


(A  -9) 
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Then  Eq.  A-8  may  be  written, 


(A-10) 


Equations  A-7  and  -10  are  the  basic  differential  equations  governing  the 
tenperature  of  the  solid  and  the  fluid  in  the  drum  type  matrix. 


A-2.6 


The  Differential  Equations  in  Dimensionless  Space  Coordinates 


New  coordinates  are  introduced: 

_2  2 

*  hS"'  „  „  “  ri 

5  •  2*K  •— s — 

p 


(A-ll) 


2  2 

..  _hS",0  2irKp0  r  ”  riN 
1  “  C"'  2  ' 


(A-12) 


where  subscript  i  refers  to  the  inside  radius  of  the  drum.  Then  con¬ 
sidering  t  =  t(  £  ,  <^)  and  ©  =  ©(  £  ,  r^),  one  can  reduce  Eq.  A-7  and  -10 
to  the  simple  forms  obtained  in  discussing  the  flat  matrix,  naxiely. 


|£  =  ©  -  t 


(A-13) 


S-— 


(A-U) 


A— 3  Reduced  Length  and  Reduced  Period 


The  reduced  length  of  the  drum  type  matrix  is 


A  =  2irK 


r2  r2 

hSm  ro  “  rl 


(A-15) 


where  tq  is  the  outer  radius  of  the  drum.  The  reduced  period  is 

„  hS'«  2wK 


(A-16) 
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1  l 


* 


The  boundary  conditions  esqpressed  in  terms  of  A  and  K  are  the 
same  for  the  drum  type  regenerator  as  for  the  flat  matrix.  Therefore, 
the  solutions  of  the  problem  of  the  flat  regenerator  becomes  directly 
applicable  to  the  solution  of  the  problem  of  the  drum  type  regenerator. 


A-h  Interpretation  of  the  Modifications  Paused  by  Curvature 

If  t  denotes  the  time  that  a  point  in  the  matrix  takes  to  move  an 
angle  0  from  the  instant  it  enters  the  fluid,  then  0  -  u>t  and,  em¬ 
ploying  Eq.  A-9,  the  definition  of  ry  becomes 


„  hs">  r  i  (T  ri," 

*1=0*-  T  “  2  (u  “ 

I-  i-i 


(A-17) 


•where  u  and  u.  are  the  radial  velocities 
i  1  r  r  • 

It  can  be  shown  that  ■£■(—  -  is  the  time 

the  fluid  to  move  radially  from  r^  to  r. 
defined  by  Eq.  A-12  has  the  same  form  as  the 
fined  by  Eq.  2-11. 


at  r  and  r^,  respectively, 
required  for  a  particle  in 
Thus  it  is  sho-wn  that 
yi  for  flat  matrices  de- 


Also,  since  2ttK/oj  represents  the  time  that  any  point  in  the  ma¬ 
trix  is  exposed  to  the  fluid,  the  expression  for  ft  in  Eq.  A-16  has 
essentially  the  same  meaning  as  the  expression  for  Tt  in  Eq.  2-19. 


Regarding  £,  ,  one  may  define  a  mean  value  of  the  mass  velocity  de¬ 
noted  by  m"  based  on  the  arithmetic  mean  of  the  facial  areas  at 
J  m,r 

radii  r.  and  r: 

l 

r  +  r. 

m'  =  m"  •  2  mr  - s— -  •  K  (A-18) 

m,r  2  v  ' 

Then 


1  =  "V  (r  -  ri)  <A-19) 

m,r  p 

Thus,  it  is  shown  that  £  defined  by  Eq.  A-ll  has  essentially  the  same 
form  as  the  £,  for  flat  matrices  defined  by  Eq.  2-10. 
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In  particular, 


a  =  rr  <ro  -  ri>  (A-20> 

in  p 

where  m”  is  the  mean  facial  mass  velocity  based  on  the  areas  at  r.  and 
m  l 

r  .  It  follows  that  to  apply  the  theory  of  the  flat  matrix  to  a  drum  re¬ 
generator,  the  reduced  length  must  be  evaluated  on  the  basis  of  the  mass 
velocity  at  the  mean  radius  of  the  drum. 


! 
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Appendix  B 

WORK  AND  THERMAL  EFFICIENCY  OF  AN  ACTUAL  REGENERATIVE  GAS  TURBINE  CYCLE 
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Appendix  Bs  WORK  AND  THERMAL  EFFICIENCY  OF  AN  ACTUAL 
REGENERATIVE  GAS  TURBINE  CYCLE 


B-l  General  Remarks 

Egressions  for  evaluating  the  ■work  and  thermal  efficiency  of  an  actual 
regenerative  gas  turbine  cycle  are  derived  in  this  appendix.  It  is  assumed 
that  the  working  medium  is  a  thermally  perfect  gas  having  constant  specific 
heats.  The  expressions  are  developed  for  a  pound  of  air  entering  the 
system  at  Poing  1  in  Fig.  4-1.  It  is  supposed  that  the  compressor  effici¬ 
ency  >|c,  the  turbine  efficiency  the  regenerator  effectiveness 

YlReg#  ^le  Pressure  ratio  =  Pj/p^j  the  leakage  fraction  Am/m,  and 
the  summation  of  the  pressure  losses  (as  well  as  their  distribution) 

2  (A p/p)  are  known. 


By  definition, 


"ST  AP  P2  -  ,  P4  -  pl 

^  P  P2  Pl 


(&-1) 


Ratios  a-^  and  are  defined: 


P4“P1 

Pl 


=  a. 


A£ 


(B-2) 


and 


P2-P3 


- - = 

Po  2 


(B-3) 


Hence, 


di  +  a2  =  1 


(B-4) 


T) ,  and 


B-2  Evaluation  of  Temperatures  T, 

It  has  been  mentioned  in  the  text  that  Points  2»  and  4r  represent 
the  ideal  end  states  of  isentropic  corpression  and  expansion,  respectively. 
Therefore, 


T2,  =  Ti  *  rp 


».  -1 
K 


(B-5) 


io4 
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V  = T 


3  r 


L  .  i-iX 


P  1  -  a. 


(-SR) 

P  ' 


(B-6) 


Tenperatures  ^  end  are  related  to  T.,,  and  T,  ,  by  the  con- 
pressor  and  turbine  efficiencies.-  ■  j  . 

1 


(B-7) 


T2  =  T1  1  + 


Ti,  =  t3  j1  +  nt  f(r  ' 1  * 

{  L\p 


K-l 


(B-8) 


Mow,  the  regenerator's  effectiveness  by  definition  is,  with  reference 
to  Fig.  U-l, 


-  TA  "  T? 

Iseg  =  -  t2 


(B-?) 


Consequently,  the  intermediate  te.per.ture  T^,  that  is,  the  teapenature 
of  the  cold  air  leaving  the  regenerator,  may  be  expressed  in  ten*  of 

T1  and  T3“ 


ic  - 1 


T  -  T  .  ^Reg  /  «■ 

1  Mv't'  p 


+  vjlifeg  1_  It-  1t^"--— 1 

\P  i-a. 


<%> 


vc-1 


(B-10) 


Having  found  Tg,  T^,  and 
now  calculate  the  heat  added 
cycle,  as  shown  below. 


m  terms  of  the  given  quantities,  one  my 
Qa  to  the  cycle  and  the  net  work  W  of  the 
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B~3  Evaluation  of  the  Heat  Added  to  the  Cycle 


The  heat  input  <2  per  pound  of  air  entering  the  system  at  Point  1 

&  A 

is  the  heat  required  to  raise  the  tenperature  of  (1  -  -~)  pound  of  the 
air  from  Point  A  to  Point  3«  Hence,  1 


«a  '  cp  t1  -  -T><T3  -  V 


(B-ll) 


B-li  Net  Work  of  the  Cycle 

The  net  work  of  the  cycle  per  pound  of  air  at  Point  1  is  the  dif¬ 
ference  between  the  work  done  by  the  turbine,  namely. 


»t  =  d  -  TT»h3  -  V  - (1  -  TT>  °p(t3  -  Tu)  .  (B-12) 

and  the  work  done  by  the  conpressor,  namely, 

»c  “  h2  -  \  =  °p(T2  -  Tl>  (B-13) 

Hence , 

W=  cp  [(1  -  -~)(T3  -  Th)  -  (T2  -  TX)J  (B-lli) 


B~5  Maximum  Thermal  Efficiency 

In  general,  the  thermal  efficiency  of  the  cycle  is 


n  =  W 

1th 


(B-1S) 


From  the  preceding  sections  of  this  appendix  it  is  evident  that  de¬ 

pends  upon  Tv  T3,  t^c,  r^.,  r^g,  2  (Ap/p),  An/m,  and  r  .  For 
any  practical  values  of  the  first  eight  quantities,  which  must  be  given 

or  calculated,  there  is  a  value  of  the  pressure  ratio  r  which  gives 

P 

the  maximum  thermal  efficiency  (cf.  Fig.  1-3).  This  value  may 

be  found  by  the  elementary  processes  of  maximizing  by  means  of  the 
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differential  calculus.  Or,  it  can  be  found  more  sinply  by  plotting 

as  a  function  of  r  . 

P 

The  maximum  work  of  the  cycle  per  pound  of  air  nay  be  found  in  similar 

■ways. 
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